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?* , ABSTRACT 

Type I X-ray bursts are thermonuclear stellar explosions driven by charged- 
particle reactions. In the regime for combined H/He-ignition, the main nuclear 
flow is dominated by the rp-process (rapid proton-captures and /3 + -decays), the 
3o;-reaction, and the ap-process (a suite of (a,p) and (p,7) reactions). The main 
flow is expected to proceed away from the valley of stability, eventually reaching 
the proton drip-line beyond A = 38. Detailed analysis of the relevant reactions 
along the main path has only been scarcely addressed, mainly in the context 
of parameterized one-zone models. In this paper, we present a detailed study 
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of the nucleosynthesis and nuclear processes powering type I X-ray bursts. The 
reported 11 bursts have been computed by means of a spherically symmetric 
(ID), Lagrangian, hydrodynamic code, linked to a nuclear reaction network that 
contains 325 isotopes (from 1 H to 107 Te), and 1392 nuclear processes. These 
evolutionary sequences, followed from the onset of accretion up to the explosion 
and expansion stages, have been performed for 2 different metallicities to explore 
the dependence between the extension of the main nuclear flow and the initial 
metal content. We carefully analyze the dominant reactions and the products 
of nucleosynthesis, together with the the physical parameters that determine 
the light curve (including recurrence times, ratios between persistent and burst 
luminosities, or the extent of the envelope expansion). Results are in qualitative 
agreement with the observed properties of some well-studied bursting sources. 
Leakage from the predicted SbSnTe-cycle cannot be discarded in some of our 
models. Production of 12 C (and implications for the mechanism that powers 
superbursts), light p- nuclei, and the amount of H left over after the bursting 
episodes will also be discussed. 

Subject headings: hydrodynamics - nuclear reactions - nucleosynthesis - abun- 
dances - stars: neutron - X-rays: bursts 



1. Introduction 

Type I X-ray bursts (hereafter, XRBs) are cataclysmic stellar events. They are powered 
by thermonuclear runaways (TNR) in the H/He-rich envelopes accreted onto neutron stars 
in close binary systems (see reviews by Bildsten 1998; Lewin et al. 1993, 1995; Psaltis 
2006; Schatz & Rehm 2006; Strohmayer & Bildsten 2006). These events constitute the most 
frequent type of thermonuclear stellar explosion in the Galaxy (the third, in terms of total 
energy output after supernovae and classical novae) , in part because of their short recurrence 
period (hours to days). About 90 Galactic low-mass X-ray binaries exhibiting such bursting 
behavior (with burst durations of T^ urst ~ 10 - 100 s) have been found since the discovery 
of XRBs by Grindlay et al. (1976), and independently, by Belian et al. (1976). Type I X- 
ray bursts and their associated nucleosynthesis have been extensively modeled by different 
groups (see pioneering work by Woosley & Taam 1976, Maraschi & Cavaliere 1977, and Joss 
1977), reflecting the astrophysical interest in determining the nuclear processes that power 
the explosion, the light curve, as well as in providing reliable estimates for the chemical 
composition of the neutron star surface (see Schatz et al. 1999, Parikh et al. 2008, and 
references therein). 
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With a neutron star hosting the explosion, temperatures and densities in the accreted 
envelope reach high values: T peak > 10 9 K, and p ~ 10 6 g cm -3 . As a result, detailed 
nucleosynthesis studies require the use of hundreds of isotopes, linked by thousands of nuclear 
interactions, extending all the way up to the SnSbTe-mass region (Schatz et al. 2001) or 
beyond (the extent of the nuclear activitjo in the XRB nucleosynthesis study of Koike et al. 
2004 reaches 126 Xe). Indeed, the extent of the rp-process in XRBs is still not clear: recent 
experimental work now shows that it will be more difficult to reach the SnSbTe-mass region 
(Elomaa et al. 2009). Because of computational constraints, XRB nucleosynthesis studies 
have been traditionally performed using limited nuclear reaction networks, truncated near 
Ni (Woosley & Weaver 1984; Taam et al. 1993; Taam, Woosley, & Lamb 1996 -all using 
a 19-isotope network), Kr (Hanawa, Sugimoto, & Hashimoto 1983 -274-isotope network; 
Koike et al. 1999 -463 nuclides), Cd (Wallace & Woosley 1984 -16- isotope network), or Y 
(Wallace & Woosley 1981 -250-isotope network). On the other hand, Schatz et al. (1999, 
2001) have carried out very detailed nucleosynthesis calculations with a network containing 
more than 600 isotopes (up to Xe, in Schatz et al. 2001), but using a one- zone approach. 
Koike et al. (2004) have also performed detailed one-zone nucleosynthesis calculations, with 
temperature and density profiles obtained from a spherically symmetric evolutionary code, 
linked to a 1270-isotope network extending up to 198 Bi. 

Until recently, it has not been possible to couple hydrodynamic stellar calculations (in 
1-D) and detailed networks. Recent efforts include Fisker et al. (2004, 2006, 2007, 2008), 
and Tan et al. (2007) (~ 300 isotopes, up to 107 Te), Jose & Moreno (2006) (2640 reactions 
and 478 isotopes, up to Te), or Woosley et al. (2004) and Heger et al. (2007) (up to 1300 
isotopes with an adaptive network). This has prompted a detailed analysis of the nuclear 
activity powering the bursts. The most detailed work to date is that of Fisker et al. (2008), 
in the context of the 1-D general relativistic hydrodynamic code AGILE (Liebendorfer et 
al. 2002), linked to a nuclear reaction network containing 304 isotopes: a thorough analysis 
of the main nuclear activity in one characteristic burst is reported (although details for a 
sequence of 5 consecutive, 'representative' bursts are also outlined). However, because of 
the specific choice of metallicity (Z = 10 -3 , for the accreted matter) and mass-accretion rate 
(M ~ 10 17 g s _1 ) adopted, the nuclear activity does not extend much beyond mass A ~ 65, 
as a result of compositional inertia effects, that quench further extension of the nuclear path. 
Hence, the flow does not reach the SnSbTe-mass region, which was suggested as a natural 



lr The nuclear activity reflects the changes in composition driven by different nuclear processes (i.e., p- 
and a-capture reactions, /3-decays, ...) that take place in the envelope at different stages of the burst. In 
this work, the extent of the nuclear activity is arbitrarily defined by the heaviest nucleus that achieves a 
mass fraction > 10 -9 . 
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endpoint in XRB nucleosynthesis studies (see Schatz et al. 1999,2001). 

Clearly, the identification of the most relevant reactions in the A ~ 65 — 100 mass region 
remains to be addressed in detail in the framework of hydrodynamic simulations. This is 
particularly relevant since, as first pointed out by Hanawa et al. (1983), proton captures on 
heavy nuclei (i.e., the reprocess) have a dramatic effect on the shape of XRB light curves. 
To this end, a new set of type I X-ray bursts have been computed with SHIVA, a 1-D, spheri- 
cally symmetric, hydrodynamic, implicit, Lagrangian code, used extensively in the modeling 
of classical nova outbursts (see Jose & Hernanz 1998). The code has been linked to a fully 
updated nuclear reaction network containing 324 nuclides and 1392 nuclear processes, a sub- 
set of that used in Parikh et al. (2008), and includes the most relevant charged-particle 
induced reactions occurring between 1 H and 107 Te, as well as their corresponding reverse 
processes. It is worth noting that the size of this network is similar (though slightly larger) 
to that adopted by Fisker et al. (2008). In order to set up the reaction rate library for 
our study, we started by adopting the proton drip line predicted by Audi et al. (2003a, 
2003b). Experimental rates are available for a small subset of reactions (adopted from An- 
gulo et al. 1999, Iliadis et al. 2001, and some recent updates for selected reactions). For all 
other reactions for which experimental rates are not available, we used the rates from the 
Hauser-Feshbach codes MOST (Goriely 1998; Arnould & Goriely 2006) and NON-SMOKER 
(Rauscher & Thielemann 2000; for details see Parikh et al. 2008). Neutron captures are 
disregarded since our early test calculations revealed that they play a minor role in XRB 
nucleosynthesis. All reaction rates incorporate the effects of thermal excitations in the tar- 
get nuclei (Rauscher & Thielemann 2000). Screening factors are taken from Graboske et 
al. (1973) and DeWitt, Graboske, & Cooper (1973). For the weak interactions, /3-delayed 
nucleon emission and laboratory decay rates (Audi et al. 2003a) have been adopted. For 
a discussion of employing stellar versus laboratory decay rates, see Woosley et al. (2004). 
It is worth noting, however, that many computed stellar decay rates ( Fuller et al. 1982a, 
1982b; Langanke & Martinez-Pinedo 2000) do not converge to their laboratory values at 
lower temperatures and densities, calling into question the model used for these calculations. 
Studies employing properly converging stellar decay rates for all isotopes relevant to XRB 
nucleosynthesis have not been performed by any group yet, and would certainly be interest- 
ing, although the results presented in this work would not be dramatically affected by their 
inclusion. 

SHIVA uses a time-dependent formalism for convective transport whenever the charac- 
teristic convective timescale becomes larger that the integration time step. Partial mixing 
between adjacent convective shells is treated by means of a diffusion equation (Prialnik, 
Shara, & Shaviv 1979). No additional semiconvection or thermohaline mixing is consid- 
ered. Models make use of Iben's (1975) opacity fits, better suited than the OPAL opacities 
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for astrophysical environments that exhibit strong variations in metallicity, as in XRB nu- 
cleosynthesis. However, plans to incorporate these more realistic opacities are currently 
underway. The adopted equation of state includes contributions from the electron gas (with 
different degrees of degeneracy; Blinnikov et al. 1996), a multicomponent ion plasma, and 
radiation; Coulomb corrections to the electronic pressure are also taken into account. 

Accretion is computed by redistributing material through a constant number of envelope 
shells (see Kutter & Sparks 1980, for details). To handle this, a tiny envelope, containing 
1.1 x 10 18 g of material (less than 1 permil of the total envelope mass accreted during the 
first bursting episode), distributed through all the envelope shells, is put initially in place(the 
influence of the number of envelope shells on burst properties will be discussed in Section 
3). The model is then relaxed using a few, very large timesteps, to guarantee hydrostatic 
equilibrium. The temperature at the bottom of the envelope barely reaches 2.7 x 10 7 K, 
whereas the density is just 1.4 x 10 3 g cm -3 (corresponding to a pressure of 5.7 x 10 18 dyn 
cm -2 ). Mass accretion and nuclear reactions are then initiated. 

Special emphasis is placed on the effect of the initial metallicity of the accreted matter 
on the main nuclear path, which in turn, will affect the final post-burst envelope composition 
and the shape of the light curves. 

The structure of the manuscript is as follows: in Section 2, we analyze the main features 
(nuclear path, nucleosynthesis, light curves, etc) of a series of 4 bursts computed in a model 
with solar-like accreted material. The effect of the resolution adopted in this model is 
discussed in Section 3. A detailed analysis of the impact of the metallicity of the accreted 
material on burst properties is given in Section 4. Finally, a comparison with previous work, 
together with a thorough analysis of the corrections posed by general relativity, are discussed 
in Section 5. 



2. Model 1 

We summarize the gross properties of a series of thermonuclear bursts driven by mass 
accretion onto a 1.4 M neutron star {Lini = 1-6 x 10 34 erg s _1 = 4.14 L ), at a rate 
M acc = 1.75 x 10~ 9 M yr _1 (corresponding to 0.08 Mem)- The composition of the accreted 
material (see Table 1) is assumed to be solar-like (X=0.7048, Y=0.2752, Z=0.02). All metals 
are initially assumed to be in the form of 14 N, following the rapid rearrangement of CNO 
isotopes that naturally occurs early in the burst (see Woosley et al. 2004). This model is 
qualitatively similar to model ZM, computed by Woosley et al. (2004) in the framework 
of the 1-D, hydrodynamic, implicit code KEPLER. This choice is made intentionally to 
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compare with previous hydro dynamic studies. Notice, however, that Woosley et al. assume 
a value of 10 km for the neutron star radius. In contrast, our model yields a value of 13.1 
km, following the integration of the neutron star structurcl from the core to its surface, in 
hydrostatic equilibrium. Differences in the neutron star size (and in turn, in surface gravity) 
may effect the strength of the explosion (mass accreted, peak temperature, nucleosynthesis, 
etc). 



2.1. First burst 

The piling up of solar-like material on top of the neutron star during the accretion stage 
progressively compresses and heats the envelope (consisting of 60 shells). Indeed, only 145 
seconds since the beginning of accretion, the temperature at the base of the envelope reaches 
Tbase = 5 x 10 7 K (with p base exceeding 10 4 g cur 3 ). 

The early nuclear activity is fully dominated by H-burning through hot CNO-cycle 
reactions, initiated by proton-captures on 14 N nuclei. At this stage (t=2327 s), the enve- 
lope achieves T& ase ~ 10 8 K {pbase ~ 6.5 x 10 4 g cm" 3 ), with an energy generation rate 
of e n «c ~ 1.2 x 10 14 erg g" 1 s _1 . The main reaction path (see Fig. [I]) is led by 15 N(p, 
a) 12 C, which powers 12 C(p, 7) 13 N(p, 7 ) 14 0(/3+) 14 N(p, 7 ) 15 0(/3+) 15 N. This suite of nuclear 
processes competes with 13 N(/3 + ) 13 C(p, 7) 14 N, and to a lesser extent, with 15 N(p, 7) 16 0(p, 
7) 17 F(/3+) 17 0(p, a) 14 N. Besides H (X=0.689) and 4 He (Y=0.290), the next most abundant 
species in the envelope is now 15 (1.2 x 10 -2 ). The amount of unburned 14 N has dropped 
to 2.4 x 10~ 3 (in the following, when discussing the nucleosynthesis, we will refer to abun- 
dances by mass, i.e., mass fractions). Other CNO-group nuclei, such as 12 C (5.4 x 10~ 5 ), 
13 N (4.2 x 10" 4 ), 14 (6.1 x 10~ 3 ), or 16 (1.1 x 10" 5 ), have already achieved an abundance 
> 10~ 5 . 

4.49 hours (16,163 s) after the beginning of accretion, Tbase reaches 2.1 x 10 8 K. Mass 
accretion in highly degenerate conditions has compressed the envelope base to a density of 
Pbase = 2.7 x 10 5 g cm -3 (Phase = 9.1 x 10 21 dyn cm" 2 ). The total luminosity of the star 
has now increased to a value of 2.5 x 10 35 erg s _1 . The main nuclear activitjo (Fig. [2]) is 



2 Thc Harrison- Wheeler equation of state (see Shapiro & Teukolsky 1983) has been adopted to build up 
the initial neutron star configuration in hydrostatic equilibrium conditions. Although this equation of state 
is a rather crude approximation to the physical conditions in the very deep interior of a neutron star, the 
radius obtained is in good agreement with the values derived with more accurate equations of state (Lattimer 
2009). 

3 Equilibrium (p,7)-(7),p) pairs have been removed from the accompanying plots to highlight those reac- 
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still dominated by proton-captures and /3 + -decays, characteristics of the hot CNO-cycle, now 
supplemented by 15 N(p, 7) 16 0(p, 7) 17 F(p, 7) 18 Ne(/3 + ) 18 F(p, a) 15 0, and by the 3a reaction. 
The numerous p-captures have reduced the hydrogen content to 0.408. In turn, 4 He has 
increased to 0.570, becoming now the most abundant species at the base of the envelope 
(followed by the short-lived species 14 [7.7 x 10~ 3 ] and 15 [1.4 x 10~ 2 ]), while most of the 
CNO nuclei have been reduced to 10 -7 — 10 -9 , by mass). The extension of the main nuclear 
activity reaches 40 Ca. Indeed, 32 S and 40 Ca are the only species in the Ne-Ca mass region 
with abundances exceeding 10~ 9 . 

Convection sets in erratically, at ~ 1 m above the core-envelope interface (the overall 
envelope size, Az, is ~ 14 m), when Tbase reaches 3.9 x 10 8 K, and progressively extends 
throughout the whole envelope. Time evolution of density, temperature, pressure, and rate 
of nuclear energy generation, at the innermost envelope shell, as well as of the overall neutron 
star luminosity and envelope size, are shown in Figs. [7]and|8j 

Shortly after, at t=5.88 hours (21,181 s), T hase reaches 4 x 10 8 K (with p base = 2.9 x 10 5 
g cm -3 , and Pbase — 1-2 x 10 22 dyn cm -2 ). The hydrogen content has dropped to 0.209, 
whereas 4 He achieves 0.625. In turn, the rate of nuclear energy generation has increased 
to a value of 2.8 x 10 16 erg g -1 s -1 . The metallicity of this innermost envelope shell has 
increased from an initial value of 0.02 to 0.17, due to leakage from CNO cycle (mainly 
powered by 15 0(a, 7)). As before, the next most abundant species are 14 (6.9 x 10~ 2 ), 
and 15 (6.5 x 10~ 2 ), but the number of isotopes with moderately large abundances has now 
increased. Indeed, 40 Ca, 22 Mg, 18 Ne, 34 Ar, 48 Cr, and 42 > 44 Ti, have achieved mass fractions of 
the order of 10 -3 . The nuclear activity extends as far as 53 Co. The largest reaction fluxes 
(number of reactions per unit time and volume) correspond to the equilibrium processes 
21 Mg(p, 7) 22 A1( 7 , p) 21 Mg, 30 S(p, 7) 31 C1( 7 , p) 30 S, and 25 Si(p, 7) 26 P(7, p) 25 Si. Additional 
activity is powered by 3a ^ 12 C(p, 7) 13 N(p, 7 ) 14 0, followed by 14 0(a, p) 17 F(p, 7 ) 18 Ne. The 
suite of secondary nuclear paths is rich and complex (see Fig. [3]) , and is mainly dominated by 
p-capture reactions and /3 + -decays, as well as by the CNO-breakout reaction 15 0(a, 7) 19 Ne. 
It is worth noting that the main nuclear path above Ca begins to move away from the valley 
of stability, towards the proton-drip line (see Fig. [3j lower panel). 

Just 2.3 seconds later (t = 21,183 s), T fease achieves 5 x 10 8 K. p base has slightly decreased 
to 2.3 x 10 5 g cm -3 because of a mild envelope expansion (Az ~ 15.5 m). Notice, however, 
that Pbase = 1-2 x 10 22 dyn cm -2 . Hence, the TNR is taking place nearly at constant 
pressure. A time-dependent, convective mixing with adjacent shells, with a characteristic 
timescale of r conv ~ 10~ 4 s (v conv ~ 10 3 — 10 5 cm s _1 ), causes a slight increase in the H 



tions of lower flux that directly lead to the production of heavier species during the burst. 



abundance at the base of the envelope. Indeed, the H abundance is now 0.288, by mass, 
whereas the 4 He content has slightly decreased to 0.563 (due to the high temperatures, 
favoring a-captures). The next most abundant species is now 18 Ne (4.4 x 10~ 2 ), together 
with 14 - 15 (4.2 x 10~ 2 , and 1.8 x 10~ 2 , respectively). Several isotopes, such as 21 - 22 Mg, 29 ' 30 S, 
50,52p e ^ 27p^ 24,25g^ 49,50,5i]y[ n ^ anc } 34^.^ j^yg achieved abundances an order of magnitude 

lower (~ 10~ 3 ). The nuclear activity extends up to 57 Cu now, powering an energy generation 
rate of 1.2 x 10 17 erg g _1 s _1 , and an overall luminosity of 1.3 x 10 36 erg s -1 . The largest 
fluxes still correspond to the forward and reverse reactions 30 S(p, 7) 31 C1(7, p) 30 S, 21 Mg(p, 
7 ) 22 A1( 7 , p) 21 Mg, and 25 Si(p, 7) 26 P(7, p) 25 Si, together with 14 0(a, p) 17 F(p, 7 ) 18 Ne. 14 0+a 
becomes the most important a-capture reaction, overcoming 15 0(a, 7), or the 3a. Additional 
activity is driven by 18 Ne(/3+) 18 F(p, a) 15 0, 19 Ne(p, 7) 20 Na(p, 7) 21 Mg, 21 Na(p, 7) 22 Mg(p, 
7 ) 23 A1, 12 C(p, 7 ) 13 N(p, 7 ) 14 0, and 26 Si(p, 7 ) 27 P. 

A qualitatively similar picture is found when Tb ase achieves 7 x 10 8 K (t = 21,185 
s), with the most abundant species at the envelope base being H (0.308), 4 He (0.507), 
18 Ne (4.8 x 10- 2 ), 22 Mg (2.5 x 10~ 2 ), 29 ' 30 S (1.1 x 10~ 2 , and 2.2 x KT 2 , respectively), and 
24,25gj ^ ^ x iq-2 anc [ 17 x io -2 , respectively). The number of species with abundances 
of the order of 10" 3 includes now 54 > 55 > 56 Ni, 15 0, 28 S, 52 Fe, 27 P, 38 Ca, and 33 - 34 Ar, with the 
main nuclear activity (see Fig. H]) extending all the way up to 60 Zn. The largest reaction 
fluxes are achieved by the equilibrium processes described before, supplemented now by 
26 Si(p, 7) 27 P( 7 , p) 26 Si, 22 Mg(p, 7) 23 A1( 7 , p) 22 Mg, 29 S(p, 7) 30 C1(7, p) 29 S, and 16 0(p, 7 ) 17 F( 7 , 
p) 16 0, followed by p-capture reactions (and /3 + -decays) on Ne-Mg nuclei, such as 19 Ne(p, 
7) 20 Na(p, 7) 21 Mg(/3+) 21 Na(p, 7) 22 Mg, or 17 F(p, 7 ) 18 Ne(/3 + ) 18 F(p, a) 15 0, and by the a- 
capture reactions 15 0(a, 7) 19 Ne, 14 0(a, p) 17 F, and the 3a. 

One second later (t = 21,186 s), T fease achieves 9 x 10 8 K. The hectic nuclear activity, 
which at this stage releases e nuc ~ 3.7 x 10 17 erg g _1 s -1 , has reduced the H and 4 He 
abundances down to 0.262 and 0.457, respectively. The next most abundant species are now 
22 Mg, 25 Si, 28 ' 29 > 30 S, 33 - 34 Ar, and 60 Zn, all with mass fractions ~ 10 -2 . The main nuclear 
activity has extended up to 68 Se. Aside from equilibrium (p, 7) and (7, p) pairs (that 
involve 16 0- 17 F and a handful of species in the mass range Mg-Zn), the largest reaction 
fluxes correspond to a suite of p-captures and /3 + -decays (see Fig. \5\ lower panel), mainly 
25 Al(p, 7) 26 Si, and 27 > 28 < 29 p(p 5 ry)28,29,30g_ Moreover, the most important a-capture reactions, 
22 Mg(a, p) 25 Al, the 3a, 18 Ne(a, p) 21 Na, and 14 0(a, p) 17 F, have fluxes of the order of log 
F ~ -3 (notice the moderate extension of a-captures towards heavier species as a result of 
the higher temperatures). A very limited nuclear activity in the A=65-100 mass region is, 
at this stage, driven by 65 Ge(p, 7) 66 As(p, 7) 67 Se(/3+) 67 As(p, 7) 68 Se (with log F ~ -8), and 
66 As(/3+) 66 Ge(p, 7 ) 67 As (log F ~ -9). 
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At t = 21,188 s, when Tbase achieves 1 x 10 9 K, the energy generation rate by nuclear 
reactions reaches its maximum value: e nuC)max ~ 4.1 x 10 17 erg g _1 s _1 . Two seconds later, 
the envelope will attain maximum expansion, with a size Az max ~ 44 m. 

Proton and a-captures continue to reduce the overall H and He abundances at the enve- 
lope base (0.191 and 0.400, respectively). The next most abundant species is now 30 S (0.103) 
-a waiting point for the main nuclear path-, followed by 33 > 34 Ar, 37,38 Ca, 42 Ti, 46 Cr, 50 Fe, 56 Ni, 
and 60 Zn (all with Xj ~ 10~ 2 ). The nuclear activity has reached 76 Sr (see Fig. E]). In the 
A=65-100 mass region, in particular, the nuclear activity is now dominated by the chains 
65 Ge(p, 7 ) 66 As(p, 7 ) 67 Se(/3+) 67 As(p, 7 ) 68 Se (with log F ~ -7), 66 As(/3+) 66 Ge(p, 7 ) 67 As (log 
F ~ -8), 65 As(p, 7 ) 66 Se(/3+) 66 As, and 68 Se(/3+) 68 As(p, 7 ) 69 Se(p, 7 ) 70 Br(p, 7 ) 71 Kr(/3+) 71 Br(p, 
7 ) 72 Kr (log F ~ -9). In terms of energy production, the most important contributions come 
from 39 Ca(p, 7 ) 40 Sc (log F ~ -2), multiple (p, 7 )-reactions involving species in the mass 
range A ~ 20-60, and a handful of a-capture reactions, such as 22 Mg(a, p) 25 Al, 18 Ne(a, 
p) 21 Na, the 3a, and 14 0(a, p) 17 F (log F ~ -3). 

Four seconds later, the envelope base achieves a maximum temperature of T pea k ~ 
1.06 x 10 9 K (similar values are reported in the simulations by Fisker et al. 2008). Besides 
H (0.220) and 4 He (0.370), the next most abundant isotope is now 60 Zn (0.159) -another 
waiting point for the nuclear flow-, followed by 30 S (3.3 x 10~ 2 ), 34 Ar (3.2 x 10~ 2 ), 38 Ca 
(2.4 x 10~ 2 ), 46 Cr (1.8 x 10" 2 ), 50 Fe (1.3 x HT 2 ), 55 ' 56 Ni (1.8 x 10" 2 , and 2.3 x 10~ 2 , re- 
spectively), and 59 Zn (1.8 x 10 -2 ). As shown in Fig. [91 the extension of the main nuclear 
path reaches 80 Zr. The nuclear activity in the A=65-100 mass region is now dominated 
by 65 Ge(p, 7 ) 66 As(p, 7 ) 67 Se(/3+) 67 As(p, 7 ) 68 Se (log F ~ -6), 65 As(p, 7 ) 66 Se(/3+) 66 As (log F 
~ -7), and 66 As(/3+) 66 Ge(p, 7 ) 67 As, 68 Se(/3+) 68 As(p, 7 ) 69 Se(p, 7 ) 70 Br(p, 7 ) 71 Kr(/3+) 71 Br(p, 
7 ) 72 Kr (log F ~ -8). Energy production is due to dozens of (p, 7 )-reactions involving species 
in the mass range A ~ 20-60, plus some a-capture reactions, such as the 3a, 14 0(a, p) 17 F, 
22 Mg(a, p) 25 Al, and 18 Ne(a, p) 21 Na (log F ~ -3). Less than a second later (t = 21,192.3 s), 
the neutron star reaches maximum luminosity, L max = 3.8 x 10 38 erg s _1 (9.8 x 10 4 L Q ). 

The numerous proton-captures on many species during the decline from T pea k reduce 
dramatically the H content in the innermost shell. Indeed, when Tf, ase achieves 9.3 x 10 8 K 
(t = 21,200 s), the H abundance drops below 0.1, while X( 4 He) = 0.283. Actually, the most 
abundant species in this shell is now 60 Zn (0.43 by mass), followed by 30 S (2.7 x 10~ 2 ), 34 Ar 
(2.1 x 10~ 2 ), 38 Ca (1.3 x 10~ 2 ), 56 Ni (1.8 x 10" 2 ), and 64 Ge (3.3 x 10~ 2 ). The nuclear activity 
reaches 90 Ru. 

Five seconds later (t = 21,205 s), when T fease drops to 9.0 x 10 8 K, 60 Zn achieves a 
maximum abundance of 0.519, by mass. H has been reduced to 1.3 x 10~ 2 (X( 4 He)=0.226). 
The next most abundant species are now 26 Si, 30 S, 34 Ar, 38 Ca, 56 Ni, and 64 Ge (all with 
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mass fractions ~ 10 -2 ). The nuclear activity has not progressed beyond 90 Ru. The largest 
reaction fluxes correspond to proton-captures and reverse photodisintegration reactions at 
equilibrium. Many other nuclear processes (/3 + -decays and a-induced reactions like the 
triple-a, 14 0(a, p) 17 F, or 22 Mg(a, p)) contribute to the overall nuclear activity (see Fig. [HIj) . 
In the A=65-100 mass region, this is driven by 65 Ge(p, 7) 66 As (followed either by 66 As(p, 
7 ) 67 Se(/3+) 67 As, or by 66 As(/3+) 66 Ge(p, 7 ) 67 As), 67 As(p, 7 ) 68 Se(/3+) 68 As(p, 7 ) 69 Se (log F 
~ -5), 69 Se(p, 7 ) 70 Br(/3+) 70 Se(p, 7 ) 71 Br(/3+) 71 Se, 74 Rb(/3+) 74 Kr (log F ~ -6), and 70 Br(p, 
7 ) 71 Kr(/3+) 71 Br, 72 Kr(/3+) 72 Br(p, 7 ) 73 Kr(p, 7 ) 74 Rb, 74 Kr(p, 7 ) 75 Rb(p, 7 ) 76 Sr (log F ~ -7). 
Energy production is due to (p, 7 )-reactions involving species in the mass range A ~ 20-55, 
and also to two a-capture reactions, 14 0(a, p) 17 F, and 22 Mg(a, p) 25 Al (log F ~ -3). 

Following the fast decline in temperature, when T base = 8.0 x 10 8 K (t = 21,212 s), 
the 60 Zn abundance has dropped to 0.509, due to /3 + -decays. H has been heavily depleted 
(5 x 10~ n ), whereas 4 He has been slightly reduced to an abundance of 0.190. The next most 
abundant species are 12 C, 30 P, 39 K, 56 Ni, 60 Cu, and 64 Ge. The extent of the nuclear activity 
(Fig. [TT]) is still restricted to 90 Ru; it will not proceed beyond this endpoint, first because 
of the heavy H depletion, and second, because the temperature is already too low to allow 
proton- or a-captures on heavier species due to their large Coulomb barriers. At this stage, 
the single, most important reaction, in terms of reaction fluxes, is the triple-a, followed by 
a suite of /3+-decay reactions, such as 26 Si(/3+) 26m Al(/3+) 26 Mg, 34 Cl(/3+) 34 S, 60 Zn(/3+) 60 Cu, 
or 27 Si(/3 + ) 27 Al. Several a-captures follow the triple-a reaction as a chain: 12 C(a, 7 ) 16 0(a, 
7 ) 20 Ne(a, 7 ) 24 Mg(a, 7 ) 28 Si(a, 7 ) 32 S, or through alternative paths, proceeding close to the 
valley of stability, up to ~ Ar, such as 13 N(a, p) 16 0, 25 < 27 Al(a, p) 28 ' 30 Si, 22 Mg(a, p) 25 Al, 
22 Na(a, p) 25 Mg, or 26,27 Si(a, p) 29 ' 30 P, to quote some representative cases (see Fig. [TT]) . e nuc 
has already declined to a value of ~ 3.8 x 10 15 erg g _1 s~ 4 . 

When Tbase drops to 4.3 x 10 8 K (t = 21,254 s), the nuclear energy generation rate has 
declined to a value of e nuc ~3x 10 14 erg g~ 4 s _1 . As shown in Fig. [T21 the nuclear activity 
is dominated by 60 Zn(/3 + ) 60 Cu, because of its very large abundance (X( 60 Zn)=0.416), and 
is followed by the triple-a reaction (X( 4 He)=0.137), and by a suite of /3 + -decays of very 
abundant isotopes, such as 64 Ge, 30 P, 64 Ga, or 60 Cu (all with mass fractions ~ 10 -2 , except 
60 Cu [0.104]), followed by those of 38 K, 26m Al, 68 Se, 25 A1, 68 As, 63 Ga, 59 Cu, and 61 Zn. Other 
species, such as 12 C, 26 Mg, 34 S, 39 K, or 56 Ni, have achieved an abundance of 10~ 2 by mass 
at this stage. The envelope has already shrunk to a size Az ~ 13 m, whereas the overall 
luminosity of the star has decreased to L^s = 7.7 x 10 36 erg s~ 4 . 

When Tbase reaches 2 x 10 8 K (t = 21,618 s), 60 Zn has remarkably decayed into 60 Cu, 
which now constitutes the most abundant species (with 0.393) at the base of the envelope. 
Because of the relatively low temperatures, the 4 He abundance is kept constant, at about 
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0.136. The next most abundant species are now 12 C, 26 Mg, 30 Si, 34 S, 39 K, 56 Ni, 60 Ni, 64 Ga, 
and 64 Zn. With the exception of the limited contribution of the triple-a reaction, the main 
nuclear path (Fig. [T3]) is fully dominated by a suite of /3 + -decays on numerous species, all 
the way up to 72 Br. 

When t = 28,250 s, T(, ase reaches a minimum value of 1.67 x 10 8 K, which we consider 
to mark the end of the first burst in our simulations. 60 Cu has decayed already into 60 Ni, 
now the most abundant species at the envelope base with a mass fraction of 0.504, followed 
by 4 He (0.136, not fully consumed during the TNR), and by 12 C, 26 Mg, 30 Si, 34 S, 39 K, 56 Ni, 
60 Cu, and 64 Zn (see Fig. Hi]) . The marginal nuclear activity played by a suite of /5-decays 
powers a rate of nuclear energy generation of e nuc ~ 8.8 x 10 11 erg g _1 s _1 . At this stage, 
the size of the envelope has shrunk to Az ~8m, whereas the overall luminosity of the star 
has decreased to L NS = 8.3 x 10 34 erg s _1 . 

Profiles of density, temperature, rate of energy generation, pressure, and size, along the 
accreted envelope, for different snapshots during the first bursting episode computed in this 
model, are shown in Figs. [T5] and [U5J 

It is worth noting that the main nucleosynthetic activity takes place at the innermost, 
hottest envelope shell. Because of their lower temperatures and densities, all layers above 
the ignition shell show a similar nucleosynthetic pattern but somewhat diluted, limiting the 
extent of the nuclear activity to lower masses. Even though the specific reaction sequences 
have a clear dependence with depth (see Fisker et al. 2008, for the corresponding analysis 
in one bursting episode), it is clear that the identification of the main nuclear processes 
responsible for the nucleosynthesis in XRBs can rely on an accurate account of the activity 
at the ignition shell, at different stages of the TNR, as performed in this paper. However, 
we would like to outline schematically how depth influences the extent of the nuclear ac- 
tivity throughout the envelope: while in the innermost shells of our computational domain 
(encompassing 6.8 x 10 20 g) the nuclear activity reaches 90 Ru, at 1.7 x 10 21 g above the 
core-envelope interface, the activity stops around 78 Sr (the final mass fraction of 90 Ru barely 
reaches 10~ 15 , by mass), while close to the surface (2.3 x 10 21 g), the nuclear activity does 
not extend beyond 72 Se (X( 90 Ru) ~ 10~ 18 ). Moreover, the purely nucleosynthetic imprint 
in these shells is difficult to assess since it is partially poisoned by changes in the chemical 
composition driven by convective mixing. 

All in all, the mean, mass-averaged chemical composition of the envelope at the end 
of this first burst, is mainly dominated by the presence of intermediate-mass elements (far 
below the SnSbTe-mass region). This includes 60 Ni (0.32), 4 He (0.31), 4 H (0.17), 64 Zn (0.03), 
12 C (0.02), 52 Fe (0.02), or 56 Fe (0.02) (see Table 2, for the mean composition of all species 
-stable or with a half-life > 1 hr- which achieve Xj > 10 -9 ), with a nucleosynthesis endpoint 
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(defined by the heaviest isotope with Xj > 10~ 9 ) around 89 Nb (in agreement with the results 
reported by Fisker et al. 2008). In terms of overproduction factors, / (ratio of the mass- 
averaged composition of a given isotope over its solar abundance; see also Figs. [T7] and 
□ED, 43 Ca, 45 Sc, 49 Ti, 51 V, 60 - 61 Ni, 63 ' 65 Cu, 64 > 67 > 68 Zn, 69 Ga, 74 Se, and 78 Kr achieve a value of 

/ ~ io 4 . 

It is important to stress that the presence of unburned H and 4 He in the envelope, at the 
end of the first burst, will have consequences for the subsequent eruptions (see Section 1.2). 
Notice, however (Fig. H7|) . that since the innermost envelope is devoid of H, the next burst 
will likely initiate well above the core-envelope interface. Moreover, the presence of unburned 
12 C, particularly in the inner envelope layers, has important implications for studies of the 
physical mechanism that powers superbursts (see Sect. 2.2). 



2.2. Second, third, and fourth bursts 

For conciseness, we will focus here on the main differences between the first and suc- 
cessive bursts computed for model 1. A first, remarkable difference is due to the so-called 
compositional inertia (Taam 1980; Woosley et al. 2004), which accounts for differences in 
the gross properties of the bursts driven by changes in the chemical content of the envelope. 
Indeed, after the first burst, the accreted matter will pile up on top of a metal-enriched 
envelope (the initial metallicity, Zj n ; ~ 0.02, has risen to a mass-averaged value of Z ~ 0.52, 
at the end of the first burst) that is devoid of H at its innermost layers. This will cause a 
shift in the location of the ignition region, progressively moving away from the core-envelope 
interface in successive flashes (see Fig. [TSJ right panel). 

This is schematically shown as well in Fig. [18] (left panel), which depicts the mass 
above the neutron star core, as well as the extension of the convective regions throughout the 
envelope, for the four bursting episodes computed in this model: the peaks of the explosions 
(~ 10 — 100 s) correspond to the flat regions of the diagram, whereas the stages of steady 
accretion (~ 5 hr) are indicated by the steep slopes. Notice that, in agreement with previous 
work (Woosley et al. 2004, Fisker et al. 2008), convection mainly develops around the peak 
of the bursts (during most of the explosion, energy transport is carried by radiation only). 
Notice also that because of fuel consumption (H, He), the location of the ignition shell (and 
the extent of the convective regions) moves progressively away in mass from the neutron star 
core. 

From the nucleosynthesis viewpoint (Table 2 and Figs. [T71 and IT9|) . the nuclear activity 
extends progressively towards heavier species, reaching endpoints (Xj > 10~ 9 ) around 89 Nb 
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(1 st burst), 97 Ru (2 nd burst), "Rh (3 rd burst), and 100 Pd (4 th burst). The main nuclear 
activity and the dominant reaction fluxes at peak temperature, for the 4 th burst computed 
in this model, are shown in Fig. [5U The overall mean metallicity of the envelope at the 
end of each burst is 0.52 (1 st burst), 0.71 (2 nd burst), 0.80 {3 rd burst), and 0.86 {4 th burst). 
This increase in Z reflects both the nuclear activity during each individual burst and the 
accumulated ashes from previous bursts. A similar mass-averaged 12 C yield of ~ 0.02 is 
systematically obtained at the end of each of the 4 bursts computed. This is not enough 
to power a superburst, which requires X( 12 C) min > 0.1, at the envelope base (see Cumming 
& Bildsten 2001, Strohmayer & Brown 2002, Brown 2004, Cooper & Narayan 2004, 2005, 
Cumming 2005, or Cooper et al. 2006). With respect to overproduction factors, the increase 
in nuclear activity reported for successive bursts translates also into larger values, as high 
as / ~ 10 6 , for 76 Se, 78 > 80 Kr, and 84 Sr, or / ~ 10 5 , for species such as 64 ' 68 Zn, 72 ' 73 Ge, 74 ' 77 Se, 
82 Kr, 86 ' 87 Sr, 89 Y, and 94 Mo, for the 4 th bursting episode (see Figs. IT71 and H91) . 

A summary of the gross properties of the four bursts computed in this model is given 
in Table 3. Peak temperatures and luminosities amount to T pea k ~ (1.1 — 1.3) x 10 9 K, and 
Lpeak ~ (1 — 2) x 1O 5 L , respectively. Recurrence times between bursts of r rec ~ 5 — 6.5 hr 
and ratios between persistent] 4 ] and burst luminosities of a ~ 35 — 40 (except for the first 
burst) have been obtained. These values are in agreement with those inferred from some 
observed XRB sources (see Galloway et al. 2008) such as the textbook burster GS 1826-24 
[r rec = 5.74 ±0.13 hr, a = 41.7 ± 1.6], 4U 1323-62 [r rec = 5.3 hr, a = 38 ±4], or 4U 1608-52 
[ T rec — 4.14 — 7.5 hr, a = 41 — 54]. As reported by Woosley et al. (2004), there is also some 
trend towards stabilization of these values with increasing burst number. 

Fig. [201 depicts the corresponding light curves from the second to the fourth burst. A 
quite interesting feature, observed in some XRBs such as 4U 1608-52 (Penninx et al. 1989), 
4U 17+2 (Kuulkers et al. 2002), or 4U 1709-267 (Jonker et al. 2004), is the appearance of 
a double-peaked burst in Fig. [2D] (lower left panel). Double (or triple) peaked bursts can 
be classified in two categories (Watts & Maurer 2007): the first one corresponds to the so- 
called photospheric radius expansion bursts, which exhibit multi-peaked bursts in the X-ray 
band but not in the bolometric luminosity. The second type of multi-peaked events are also 
visible in the bolometric light curves, and have been attributed to different causes, such as a 
stepped release of thermonuclear energy caused either by mixing induced by hydrodynamic 
instabilities (Fujimoto et al. 1988) or driven by a nuclear waiting-point impedance in the 



4 We define a = J. r " c L(t) dtj J , ° 01 Lit) dt, with the latter term integrated over the time during which 
the burst exceeds 1% of its peak luminosity, ro.oi- Notice that during the interburst period, the accretion 
luminosity, L acc — GMM/R ~1.5x 10 37 ergs" 1 , will hide the thermal emission from the cooling ashes, as 
shown in Figs. [8land[20l 
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thermonuclear reaction flow (Fisker et al. 2004). A preliminary analysis of the 4 th burst 
reported for model 1 suggests a likely nuclear physics origin (waiting-point impedance) for 
this double-peaked feature (see Jose & Moreno 2010). 



3. Model 2 

In the previous Section, we have reported results from a sequence of type I X-ray bursts 
computed with a coarse resolution, in which the accreted envelope was discretized in 60 
shells. We have checked the influence of the adopted number of envelope shells on the 
gross properties of the bursts by performing another simulation, identical to model 1, but 
computed with a finer resolution: 200 shells (hereafter, model 2). 

A summary of the main properties of the two bursts computed for model 2 is given in 
Table 4: the recurrence times obtained are in the same range as those reported for model 1 
(4-6 hr). The same applies to the ratios between persistent and burst luminosities, as well 
as to peak temperatures and luminosities. Similar light curves have also been obtained. 

There is also good agreement from the nucleosynthesis viewpoint, with only minor 
differences in the final, mass-averaged abundances, as shown in Table 5 (particularly, for the 
heavier species of the network, since their low abundances are very sensitive to the specific 
thermal history of the explosion; see also Fig. [T9|) . It is worth noting that both models reach 
almost identical nucleosynthesis endpoints (Xj > 10~ 9 ): 89 Nb, for the first burst computed 
in both models, and 97 Ru (model 1) and "Rh (model 2), for the second burst. Furthermore, 
the amounts of unburned H, 4 He, and 12 C are very similar in both models. As expected from 
the abovementioned similarities, there is also good agreement in terms of overproduction 
factors, dominated by 64 Zn and 60 Ni (with / ~ 10 4 ) in the first burst, and by 64 Zn, 72 Ge, 
74,76 Se, 78 ' 80 Kr, and 84 Sr (/ ~ 10 5 ), in the second, for both models. 

All in all, we conclude that the resolution adopted in model 1 is appropriate for XRB 
simulations. This is in agreement with the studies performed by Fisker et al. (2004), who 
concluded that the minimum discretization of the envelope, in 1-D hydrodynamic simulations 
of X-ray bursts, is about 25 shells. 



4. Model 3 

To test the impact of the metallicity of the accreted material (which reflects the surface 
composition of the companion star) on the overall properties of the bursts, we have computed 
another series of bursts (hereafter, model 3), driven by accretion of metal-deficient material 
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(Z ~ Z Q /20) onto a 1.4 M neutron star (L ini = 1.6 x 10 34 erg s" 1 = 4.14 L ), at a rate 
M acc = 1.75 x 10 -9 M yr _1 . The composition of the accreted material is assumed to be 
X=0.759, Y=0.240, and Z=10~ 3 , and as for model 1, all metals are initially added up in the 
form of 14 N. This model is indeed qualitatively similar to model zM, from Woosley et al. 
(2004) (see also Fisker et al. 2008, and Schatz et al. 2001). Both the envelope zoning and 
the initial relaxation phase are identical to those described for model 1. 



4.1. First burst 

The piling up of matter on top of the neutron star during the accretion stage progres- 
sively compresses and heats the envelope. 

At t=4337 s, the envelope achieves Tbase ~ 10 8 K (pbase ~ 9.8 x 10 4 g cm -3 ). The 
nuclear activity is fully dominated by the CNO-cycle. In contrast to model 1, the smaller 
metallicity of this model limits substantially the role of proton captures. Indeed, at this 
stage, H has only been reduced to 0.757 at the envelope base. The main reaction fluxes 
are actually an order of magnitude lower than those reported from model 1, for the same 
temperature, powering an energy generation rate of e nuc ~ 6.2 x 10 12 erg g" 1 s -1 . Because 
of the lower metallicity of this model, the time required to achieve Tbase ~ 10 8 K is about 
twice the value reported for model 1, resulting in a thicker, more massive envelope which will 
affect the forthcoming explosion. Besides H and He (0.242), by far the most abundant nuclei 
in the envelope, the nuclear activity in the CNO region increases the chemical abundances of 
many species in this mass range, with 15 (6.3 x 10~ 4 ) being the most abundant CNO-group 
nucleus at the envelope base. 

16.8 hours (60,347 s) after the beginning of accretion, Tbase has reached 2 x 10 8 K (with 
Phase = 5.7 x 10 5 g cm -3 , and Pbase = 3.4 x 10 22 dyn cm -2 ). The total luminosity of the star 
is only 6.9 x 10 34 erg s -1 . The main nuclear activity is governed by 15 N(p, a) l2 C and other 
reactions of the CNO cycle. 7 Be(p, 7) 8 B is at equilibrium with its reverse photodisintegration 
reaction 8 B(7, p) 7 Be. Because of the limited number of CNO-catalysts in this low-metallicity 
model, some proton-proton chain reactions, such as the pep, 3 He(o;, 7) 7 Be, or 8 Be — > 2 4 He, 
are relatively frequent. In terms of chemical abundances, the now frequent p-captures have 
reduced the hydrogen content down to a value of 0.695 (while X( 4 He) = 0.303). The next 
most abundant isotopes in the network are the short-lived species 15 (1.4 x 10~ 3 ) and 14 
(8 x 10~ 4 ). The nuclear activity (Xi > 10~ 9 ) reaches 40 Ca at this stage. 

18.1 hours (65,081 s) from the onset of accretion, Tbase reaches 4 x 10 8 K. Hydrogen 
continues to decrease smoothly (X(H) = 0.689), whereas the 4 He abundance reaches 0.280. 
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The next most abundant nuclei are 14 (1.4 x 10~ 2 ), and 15 0, now followed by 52 Fe and 
18 Ne, with mass fractions of the order of 10 -3 . The nuclear activity extends up to 58 Cu. The 
largest reaction fluxes correspond to different processes that operate almost at equilibrium 
with their inverse photodisintegration reactions, such as 21 Mg(p, 7) 22 A1(7, p) 21 Mg, 30 S(p, 
7 ) 31 C1( 7 , p) 30 S, 25 Si(p, 7) 26 P( 7 , p) 25 Si, and 7 Be(p, 7 ) 8 B( 7 , p) 7 Be. 

6 seconds later (t = 65,087 s), T base achieves 5 x 10 8 K. p base has slightly decreased to 
4.3 x 10 5 g cm -3 because of a mild envelope expansion (Az ~ 19.7 m). 4 He has slightly 
decreased to 0.267, as a result of the frequent a-captures driven by the high temperatures 
achieved. The next most abundant species are now 18 Ne (10~ 2 ), together with 14 > 15 (3.5 x 
10~ 3 and 7.5 x 10~ 3 , respectively), 52 Fe (2.8 x 10~ 3 ), and 34 Ar (1.1 x 10~ 3 ). The nuclear 
activity reaches 61 Ga, powering an energy generation rate of 2.8 x 10 16 erg g _1 s _1 . The 
overall stellar luminosity is now 2.5 x 10 35 erg s _1 . 

A qualitatively similar picture is found when Tbase achieves 7 x 10 8 K (at t = 65,090 s), 
with the nuclear activity extending all the way up to 68 Se. 

One second later (t = 65,091 s), T base achieves 10 9 K. The nuclear activity (with e nuc ~ 
1.5 x 10 17 erg g _1 s _1 ; see Fig. [22]) continues to reduce the H and 4 He abundances down 
to 0.680 and 0.224, respectively. The next most abundant species are now 28 > 29 > 30 S, 33 > 34 Ar, 
25 Si, 60 Zn, and 38 Ca (all with mass fractions ~ 10 -2 ), with the main nuclear path reaching 
72 Kr. The largest absolute fluxes are achieved by nuclear interactions between equilibrium 
(p, 7)- (7, p) pairs, which do not contribute to the net energy balance. Instead, the most 
important contributors to the energy production at this stage are 25 Al(p, 7) 26 Si, 27 > 29 > 30 p(p 7 

7 )28,30,31 S; 28g i(p; 7 )29p ; 32,33 C1 ( P) ^33,34^ 32,33 C1 ( P) ^33,34^ 35 Ar ( pj 7 )36 K; 35,36,37 K ( P) 

7 ) 36 > 37 ' 38 Ca, the chain 3a -> 12 C(p, 7) 13 N(p, 7) 14 0(a, p) 17 F(p, 7 ) 18 Ne(a, p) 21 Na(p, 7) 22 Mg, 
and a suite of /3+-decays, such as 25 Si(/3+) 25 Al, 28,29,303^+ )28,2 9) 30 P) and 33 Ar(/3+) 33 Cl. The 
activity in the A=65-100 mass region is dominated by the suite of reactions depicted in 
Fig. |22l lower panel, mainly 65 Ge(p, 7 ) 66 As, 67 Se(/3+) 67 As(p, 7 ) 68 Se, 65 As(p, 7 ) 66 Se(/3 + ) 66 As 
(with log F - -7), and 65 Se(/3+) 65 As (log F ~ -8). 

In contrast to model 1, which achieved a peak temperature of 1.06 x 10 9 K, model 3 
reaches relatively higher values. Hence, at t = 65,093 s, T(, ase achieves 1.2 x 10 9 K. The H and 
4 He abundances have been reduced to 0.648 and 0.205, respectively. The next most abundant 
nucleus is still 30 S (3.1 x 10~ 2 ), followed by 38 Ca (2 x 10~ 2 ), and by a large number of species 
with abundances of the order of 10" 3 , such as 36 - 37 Ca, 28 ' 29 S, 32 . 33 - 34 Ar, 58,59,60 Zll; 62,63,64^ 

48,49,50p e; 53,54,55 Ni; 41 T ^ and 44,45,46Q r> At thig gtage; the main nuc l ear pat h reac h eS 76 Sr. 

One second later, at t = 65,094 s, the rate of nuclear energy generation achieves a 
maximum value of e nuCtmax ~ 2.1 x 10 17 erg g _1 s _1 . 
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At t = 65,095 s, while Tbase = 1.3 x 10 9 K, the main nuclear path reaches 80 Zr. Because 
of the large temperature achieved, the number of proton- and a-captures increases, which in 
turn efficiently reduces the H (0.621) and 4 He (0.197) abundances. The next most abundant 
nucleus is 60 Zn (2.8 x 10~ 2 ), followed by 64 Ge, 38 Ca, 30 S, 55 Ni, and 59 Zn (with X* ~ 10~ 2 , see 
Fig. 123]) . The nuclear activity in the A=65-100 mass region is now dominated by 65 As(p, 
7 ) 66 Se(/3+) 66 As(p, 7 ) 67 Se (log F ~ -4), 65 Ge(p, 7 ) 66 As, and 67 Se(/3+) 67 As(p, 7 ) 68 Se (log F ~ 
-5; see Fig. [23]). Energy production is mainly due to suite of (p, 7 ) reactions and /3 + -decays 
involving nuclear species in the mass range A=30-62. 

Shortly after, at t = 65,098 s, a peak temperature of T pea k = 1.4 x 10 9 K is achieved at 
the envelope base. This is followed, less than a second later, by a maximum expansion of the 
envelope, Az max ~ 73.9 m, and by a maximum luminosity, L max = 4.0 x 10 38 erg s _1 (10 5 L ). 
The main nuclear path reaches 93 Pd (already beyond the nucleosynthesis endpoint achieved 
in model 1). With respect to the chemical abundances, the envelope base is still dominated 
by H (0.560) and 4 He (0.175), with 60 Zn reaching a mass fraction of 0.111. The next most 
abundant species are 64 Ge (6.3 x 10~ 2 ), and 68 Se (2 x 10~ 2 ). Regarding the activity in the 
A>65 mass region, at this stage is dominated by 66 Se(/3 + ) 66 As(p, 7 ) 67 Se(/3 + ) 67 As(p, 7 ) 68 Se 
(log F ~ -4), 65 Ge(p, 7 ) 66 As, and 68 Se(/3+) 68 As(p, 7 ) 69 Se(p, 7 ) 70 Br(p, 7 ) 71 Kr(/3+) 71 Br(p, 
7 ) 72 Kr (log F ~ -5; see Fig. [2U for additional processes down to log F ~ -8). The most 
important contributors to the energy production are at this stage 44 ' 45 V(p, 7 ) 45 ' 46 Cr, 49 Mn(p, 
7 ) 50 Fe(/3+) 50 Mn(p, 7 ) 51 Fe, 52 Co(p, 7 ) 53 Ni, 53 Co(p, 7 ) 54 Ni(/3+) 54 Co(p, 7 ) 55 Ni(/3+) 55 Co(p, 
7 ) 56 Ni, and 58 Zn(/3+) 58 Cu(p, 7 ) 59 Zn(/3+) 59 Cu(p, 7 ) 60 Zn (log F ~ -3). 

At t=65,110 s, following the decline from peak temperature, the envelope base achieves 
Tbase = 1.3 x 10 9 K (Fig. [25]). At this stage, the main nuclear activity has already reached 
the SnSbTe-mass region ( 104 Sn, in particular). The chemical abundances at the envelope 
base are still dominated by H (0.471), now followed by 64 Ge (0.162), and 68 Se (0.161), 
while 4 He has dropped to 0.131. The next most abundant species shift to 60 Zn (2.1 x 
10~ 2 ), and 72 Kr (1.9 x 10~ 2 ), with a suite of nuclei reaching ~ 10~ 3 ( 30 S, 67 Se, 37 ' 38 Ca, 
76 Sr, 62 ' 63 Ge, 59 Zn, 55 Ni, 34 Ar, or 50 Fe). The nuclear activity in the A>65 mass region is 
now powered by 65 As(p, 7 ) 66 Se(/3+) 66 As(p, 7 ) 67 Se(/3+) 67 As(p, 7 ) 68 Se(/3+) 68 As(p, 7 ) 69 Se(p, 
7 ) 70 Br(p, 7 ) 71 Kr(/3+) 71 Br(p, 7 ) 72 Kr (log F ~ -4), 65 Ge(p, 7 ) 66 As, 69 Br(p, 7 ) 70 Kr(/3+) 70 Br, 
72 Kr(/3+) 72 Br(p, 7 ) 73 Kr(p, 7 ) 74 Rb(p, 7 ) 75 Sr(/3+) 75 Rb(p, 7 ) 76 Sr, 78 Y(p, 7 ) 79 Zr, and 86 Tc(p, 
7 ) 87 Ru (log F ~ -5; see Fig. [25]). Energy production is not due to a handful of nuclear 
processes but to dozens of different reactions (from 3a — > 12 C all the way to 71 Br(p, 7 ) 72 Kr). 

Twenty-two seconds later (t=65,132 s), the temperature at the envelope base has de- 
creased to T base = 1.2 x 10 9 K (Fig. [26]). The H content has been slightly reduced to 0.370, 
whereas 4 He reaches 8.74 x 10~ 2 . Indeed, after H, the most abundant species at the envelope 
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base are now 68 Se (0.270), and 64 Ge (0.104), followed by 72 Kr (8.7 x 10~ 2 ), 76 Sr (2.8 x lO" 2 ), 
and 80 Zr (1.2 x 10 -2 ). At this stage, some of the heaviest species of the network have already 
achieved an abundance of ~ 10~ 3 (such as 88,89 Ru, 92 > 93 Pd, 96 Cd, "In, or 101 ' 102 Sn, together 
with the lighter nuclei 30 S, 38 Ca, 59 > 60 Zn, and 67 Se), with the nuclear activity extending all 
the way up to 107 Te. The nuclear activity in the A=65-100 mass region is similar to that 
described above for Tb ase = 1.3 x 10 9 K, and is depicted in Fig. [26] lower panel. 

At t=65,264 s, the envelope base achieves T base = 10 9 K (Fig. I2"7j) . Now, the most 
abundant element at the envelope base is 105 Sn (0.228). This is followed by a large number 
of species with abundances ~ 10" 2 , such as 104 Sn, 68 Se, 72 Kr, 104 In, 94 Pd, 64 Ge, 103 In, 76 Sr 
-all more abundant than H (1.8 x 10~ 2 ) and 4 He (2.6 x 10~ 2 ), at this stage-, together with 
102 > 103 Sn, 95 Ag, 107 Te, 100,101,102!^ so^ 60 Zn? and 98,99 Cd _ Notice th& ^ since the heaviest 

element included in our network, 107 Te, achieved already an abundance of 2.2 x 10~ 2 , leakage 
from the SnSbTe-mass region cannot be discarded. Notice, however, that a-emission for 
107 Te was not included, which may account for the high abundances reported here. Further 
studies to explore possible nucleosynthesis beyond the SnSbTe-mass region are underway 
with a larger network, so detailed abundances in this region should be taken with caution. 
The set of equilibrium (p, 7)-(7, p) pairs is now accompanied by a handful of /3 + -decays (such 
as 80 Zr(/3+) 80 Y, 76 Sr(/3+) 76 Rb, 84 Mo(/3+) 84 Nb, and 82 Nb(/3+) 82 Zr) as the nuclear processes 
with largest absolute fluxes, since H depletion and the low temperature limit the extent of 
charged-particle reactions. Indeed, these weak interactions will become progressively more 
important during the last stages of the burst. At this stage, the activity in the A=65- 
100 mass region is dominated by 76 Sr(/3+) 76 Rb(p, 7) 77 Sr(p, 7 ) 78 Y, 79 Y(p, 7 ) 80 Zr(/3 + ) 80 Y(p, 
7 ) 81 Zr(p, 7 ) 82 Nb(/3+) 82 Zr(p, 7 ) 83 Nb(p, 7 ) 84 Mo(/3+) 84 Nb(p, 7 ) 85 Mo, 89 Ru(p, 7 ) 90 Rh (log F ~ 
-4), and more than 60 different nuclear processes with log F ~ -5 (see Fig. [271 lower panel), 
involving nuclei in the mass range A= 65-104. Indeed, energy production is driven by (p, 7 ) 
and /3 + processes involving species in this mass range. 

At t=65,362 s, the temperature at the envelope base has already declined to T base = 
7.6 x 10 8 K (Fig. [28]). H is now fully depleted (7.6 x 10~ 12 ), while 4 He barely reaches 
1.8 x 10~ 2 . As before, the most abundant element at the envelope base is 105 Sn (0.251), 
followed by 104 In (0.142), and by a large number of species with abundances ~ 10 -2 . The 
depletion of H dramatically reduces the fluxes of most of the (p, 7 ) reactions, which are now 
overcome by many /3 + -decays (such as those affecting 68 ' 69 Se, 68 As, 64 Ge, 71 ' 72 Br, 60 Zn, 60 Cu, 
72 ' 73 Kr, 104 Sn, 64 Ga, 82 Zr, or 76 Sr), and by a suite of a-capture reactions, such as 3a — » 12 C(a, 
7 ) 16 0(«, 7 ) 20 Ne(«, 7 ) 24 Mg(«, 7 ) 28 Si(«, 7 ) 32 S, or 13 N(«, p) 16 0. 

At t=69,715 s, after a long decline, a minimum temperature is achieved at the enve- 
lope base, Tf, ase = 2 x 10 8 K (Fig. [29]), which we consider to mark the end of the first 
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bursting episode for this model. At this stage, H is fully depleted (3.2 x 10 -23 ) at the 
envelope base, while 4 He has a mass fraction of 1.5 x 10~ 2 only. The distribution of the 
most abundant elements almost follows the one described for T^ ase = 7.6 x 10 8 K (Fig. 
|2"81 . and is dominated by 105 Sn (0.251), followed by 104 In (0.147), and by a large number of 
species with abundances - 10" 2 , such as 94 Pd, 100,101,102,103^ 68 Ge; 64 Zllj 72 Se> 76 Rb; i07 Te; 

60 Ni, "Cd, 97,98 Ag, 89 ' 90 Ru, and 80 Y. At this stage, the dominant interactions are all f3 + - 
decays ( 60 < 61 Cu(/3+) 60 > 61 Ni, 66 > 67 Ge(/3+) 66 ' 67 Ga, 65 Ga(/3+) 65 Zn, 51 Mn(/3+) 51 Cr, 52 Fe(/3+) 52 Mn, 
63 Zn(/3+) 63 Cu, 56 Ni(/3+) 56 Co, or 43 Sc(/3+) 43 Ca), except for the triple-a reaction. 

Depth also influences the extent of the nuclear activity throughout the envelope, but in 
contrast to Model 1, the nuclear activity in all shells of our computational domain essentially 
reach the SnSbTe-mass region. Indeed, the inner part of the envelope (encompassing 3.4 x 
10 21 g) is, at the end of the burst, dominated by large amounts of 105 Sn and 104 In, the most 
abundant nuclei with mass fractions ~ 0.1 - 0.2; at 5.6 x 10 21 g above the core-envelope 
interface, the most abundant isotopes are H (0.26) and 4 He (0.12), while the most abundant 
species in the SnSbTe-mass region achieve a mass fraction ~ 10~ 3 ; and close to the surface 
(7.7 x 10 21 g), shells are largely dominated by the presence of unburned H (0.75) and 4 He 
(0.24), with X(SnSbTe) ~ lO" 7 . 

The mean, mass-averaged chemical composition of the whole envelope, at the end of the 
first bursting episode, is dominated by the presence of unburned H (0.18) and 4 He (0.084), 
followed by 105 Ag (0.075), 104 Pd (0.053), 64 Zn (0.042), 95 Ru (0.031), 68 Ge (0.028), 94 Tc 
(0.026), and 103 Ag (0.026), (see Table 8, for the mean composition of species -stable or with 
a half-life > 1 hr- which achieve Xj > 10~ 9 ), with a nucleosynthesis endpoint around 107 Cd. 
In contrast, the first burst computed in model 1 yielded, in general, lighter nuclei, 60 Ni, 4 He, 
1 H, 64 Zn, 12 C, and 52 ' 56 Fe, with a more modest nucleosynthesis endpoint around 89 Nb (Table 
2). 

In terms of overproduction factors, / (Fig. [30]) . while model 1 showed moderate values 
(/ ~ 10 4 ) for a handful of intermediate-mass elements, such as 43 Ca, 45 Sc, 49 Ti, 51 V, 60,61 Ni, 
63 ' 65 Cu, 64 ' 67 > 68 Zn, 69 Ga, 74 Se, or 78 Kr, model 3 achieves moderate overproduction factors 
(> 10 4 ), for all stable species heavier than 64 Zn, and as high as ~ 10 8 for 98 Ru, 102 ' 104 Pd, 
and 106 Cd. 



4.2. Second, third, fourth, and fifth bursts 

Table 7 summarizes the most relevant properties that characterize the five bursting 
episodes computed for model 3. Recurrence times between bursts of r rec ~ 9 hr (except for 
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the first one, for which r rec ~ 18 hr), ratios between persistent and burst luminosities of 
a ~ 20 — 30, and peak luminosities around L pea k ~ 1O 5 L represent the basic observables 
associated with this model. Indeed, the recurrence times obtained are in agreement with the 
values reported for the XRB sources (see Galloway et al. 2008) 1A 1905+00 [r rec = 8.9 hr], 
4U 1254-69 [r rec = 9.2 hr], or XTE J1710-281 [r rec = 8.9 hr, a = 22 - 100], A striking result 
is the quick stabilization of the recurrence times, that show a regular periodicity after the 
second burst. 

It is worth noting that both the recurrence periods and the ratios between persistent and 
burst luminosities are larger than those reported for model 1 (see Table 6, for comparison), 
showing a clear dependence on the metallicity of the accreted material: the smaller the metal 
content, the larger the recurrence time (and the smaller the value of a). 

The corresponding light curves (see Fig. I3T|) exhibit, in turn, a clear pattern: as shown 
in Fig. [32] (left panel), where light curves of the third bursting episode computed in models 1 
& 3 are compared, explosions in metal-deficient envelopes (such as model 3) are characterized 
by lower peak luminosities and longer decline times. A similar pattern has been reported by 
Heger et al. (2007), in the framework of 1-D, hydrodynamic models of XRBs performed with 
the KEPLER code. It is worth noting that no double-peaked bursts have been obtained in 
model 3. 

Larger peak temperatures, around T pea k ~ (1.3 — 1.4) x 10 9 K, have also been obtained 
in model 3. This, together with the longer exposure times to high temperatures (driven by 
the slower decline phase) cause a dramatic extension of the main nuclear path towards the 
SnSbTe-mass region or beyond. 

From the nucleosynthesis viewpoint, and as shown in Table 8 and Figs. EH & E2J the 
nuclear activity already reaches the end of the network ( 107 Te) at the late stages of the first 
bursting episode. The overall mean metallicity of the envelope at the end of each burst is 
now 0.74 (I s4 burst), 0.85 (2 nd burst), 0.90 (3 rd burst), 0.92 (4 th burst), and 0.93 (5 th burst). 
Notice that, although the accreted material is more metal- deficient in model 3 than in model 
1, the post-burst mean metallicity of the envelope is larger in model 3. This results from the 
combination of higher temperatures and longer burst durations, which favors the extension 
of the nuclear activity: for short bursts (like those obtained in metal-rich envelopes), only 
the fastest p- and/or a-capture reactions can naturally occur (those that proceed with a 
characteristic time shorter than the overall exposure time to high temperatures); in con- 
trast, for long-duration bursts, the overall number of p- and/or a-capture reactions increases 
dramatically. This, in particular, affects CNO-breakout through 15 0(a, 7) 19 Ne and 14 0(a, 
p) 17 F, which are favored in the longer bursts obtained for model 3. It is also worth noting 
that the abundance pattern obtained after the different bursts is very similar. This is clearly 
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shown in Fig. [331 that depicts the main nuclear activity at peak and at the end of the fourth 
burst computed for Model 3, when compared with that corresponding to the first burst - 
Figs. |23]&[2nj This fact justifies our emphasis on the reaction sequences that characterize 
the first burst (see discussion in Sect. 5.2). 

A final 12 C yield of ~ 2 x 10 -3 is obtained at the end of each bursting episode (except in 
the first one, for which X( 12 C)=8 x 10 -4 ). As reported for model 1, the amount of unburned 
12 C left over turns out to be too small to power a superburst. Finally, huge overproduction 
factors (see Figs. EH & I3"2"l) . involving heavy species such as 102 . 104 ^ 105 p ( i ) 98 Ru, or 94 Mo 
(with / ~ 10 8 ) have been obtained in model 3, in contrast with the somewhat more modest 
values achieved in model 1, where maximum overproduction factors are about / ~ 10 6 , and 
involving lighter species, such as 76 Se, 78 > 80 Kr, or 84 Sr. 



5. Discussion 

5.1. General relativity corrections 

The calculations reported here have been performed assuming Newtonian gravity. Since 
the envelope layers are very thin, it is easy to introduce general relativity corrections to 
this Newtonian framework (see Ayasli & Joss 1982, Lewin et al. 1993, Taam et al. 1993, 
Cumming et al. 2002, and Woosley et al. 2004). To this end, the surface gravity is rewritten 
as g = GM*/R 2 {1 + z), where M* is the mass, R* is the stellar radius (defined in such 
a way that the surface area is 47ri? 2 ), and z is the gravitational redshift given by 1 + z = 
(l — 2GM*/R*c 2 y 1 / 2 . Our models of M* = 1.4 M require R* = 14.3 km, and a gravitational 
redshift of z = 0.19. 

Following Woosley et al. (2004), once the redshift and radius are determined, it is 
straightforward to derive the set of correcting factors to the physical magnitudes described 
above for a suitable observer at infinity. Hence, recurrence times and burst durations should 
be increased by a factor 1 + z. The mass-accretion rate as well as the burst luminosity 
have to take into account both the difference in surface area (compared to the Newtonian 
framework) and the gravitational redshift term. The energy and rest mass-accretion rate 
scale as R 2 /(l + z), while the luminosity oc R 2 /(l + z) 2 . However, when M* is taken exactly 
as Mns (Newtonian framework), the surface area and redshift corrections for energy and 
mass accretion rate cancel out, since g oc (l + 2;)/R 2 = const, and hence, no correction to the 
observed burst energy or mass-accretion rate is necessary, while the luminosity correction is 
simply given by 1/(1 + z) = 0.84. In addition, the accretion luminosity for an observer at 
infinity changes only by a factor 1.012, that is, the ratio between gravitational energy released 
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per unit mass in general relativity, c 2 .z/{l + z), and the Newtonian value, GM^s/Rns- 
Finally, the luminosity measured at infinity will be smaller by a factor of (1 + z) — 1.19. 



5.2. Comparison with previous work 

For consistency, the results discussed in this paper have been compared with those 
reported in previous work (obtained with similar hydrodynamic codes or in the framework 
of one- zone models). 

As emphasized in Section 2, model 1 is qualitatively similar to model ZM of Woosley 
et al. (2004). The twelve bursts computed by Woosley et al. (2004) in a Newtonian frame 
were characterized by recurrence times of about ~ 2.7 hr, peak luminosities of L pea k ~ 
(1.5 — 2) x 10 38 erg s _1 , and ratios between persistent and burst luminosities of a ~ 60 — 65. 
Our calculations (model 1, Newtonian frame) yield r rec ~ 5 — 6.5 hr, L pea k ~ (3 — 7) x 10 38 
erg s^ 1 , and a ~ 35 — 40. 

The role played by the metallicity of the accreted material (model 3, with Z = Z Q /20 
= 0.001) qualitatively agrees with the pattern reported by Woosley et al. (2004) (see also, 
Heger et al. 2007). Longer recurrence times of ~ 9 hours, peak temperatures of about 
(1.3 — 1.4) x 10 9 K, and ratios between persistent and burst luminosities of a ~ 20 — 30 
(with Lp ea k ~ 10 38 erg s _1 ) have been obtained in the 5 bursts computed in model 3. In 
turn, the fifteen bursts computed by Woosley et al. (2004) for model zM are characterized 
by recurrence times of about 3 — 3.5 hr, peak luminosities of L pea k ~ 10 38 erg s -1 , and ratios 
between persistent and burst luminosities of a ~ 50 — 60. Results reveal a dependence of 
burst properties on the metallicity of the accreted material: the smaller the metal content, 
the larger the recurrence time (and the smaller the a). In turn, explosions in metal-deficient 
envelopes (i.e., model 3) are characterized by lower peak luminosities and longer decline 
times, in agreement with the pattern described in Woosley et al. (2004) and Heger et al. 
(2007). Model 3 bears as well a clear resemblance with the model computed by Fisker et al. 
(2008). In that work, five representative bursting sequences were analyzed, with r rec ~ 3.5—4 
hr, L pea k ~ (7 — 8) x 10 37 erg s -1 , and a ~ 65 — 70, as measured at infinity. 

Despite the qualitative similaries in the gross properties of the bursts presented in this 
paper (as well as in the role played by the metallicity of the accreted material) and those 
reported in previous work, a quantitative comparison reveals some discrepancies that are 
worth analyzing. In model 1 (with Z = Z Q ), our computations yield systematically larger 
(by a factor of ~ 2) recurrence times and peak luminosities (and hence, lower a) than model 
ZM of Woosley et al. (2004). Similar results are found in the low-metallicity case (model 3, 
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with Z = Z Q /20) when compared with model zM of Woosley et a!. (2004), except for the 
peak luminosities that turn out to be very similar. It is also worth noting that the values 
reported by Fisker et al. (2008) show discrepancies with respect to Woosley et al. (2004), 
in particular, lower peak luminosities (and larger a). A major difference concerns the much 
larger effect played by the metallicity of the accreted material in this work as compared with 
Woosley et al. (2004), who explained the moderate effect found as due to compositional 
inertia washing out the influence of the initial metallicity. Another striking issue concerns 
the extremely large differences in the gross physical characteristics -nucleosynthesis, energies 
or recurrence times- between the first and subsequent bursts, as reported by Woosley et al. 
(2004). In terms of nucleoynthesis or nuclear activity, Figs. [301 EH & (32J reveal a similar 
behavior for the different bursts (although a somewhat lower production of intermediate- 
mass elements as well as of the heaviest elements is reported for the first burst computed in 
model 3). 

Very limited information on the nucleosynthetic yields obtained in model ZM is given 
in Woosley et al. (2004). Thus, we will restrict the discussion on the extent of the nuclear 
activity and on the resulting chemical abundance pattern to model 3, through a brief com- 
parison with the work reported by Schatz et al. (2001), Fisker et al. (2008) and Woosley 
et al. (2004) (for model zM). It is worth noting that both the nucleosynthetic end-point 
(located in the SnSbTe-mass region) and the main nuclear path in the A ~ 50 — 100 mass 
region obtained in this work (passing through a suite of different nuclei such, as 55 Co, 60 Zn, 
70 Br, 75 Rb, 85 Mo, 90 Rh, or 100 - 105 Sn) are similar to those reported by Schatz et al. (2001) 
in the framework of one- zone calculations. Whereas the main nuclear path in the first burst 
of model zM (Woosley et al. 2004) is very similar to the one reported in this work, com- 
positional inertia causes a more limited extension of the nuclear activity in the successive 
bursts of Woosley et al.: hence, while the three most abundant nuclei at the bottom of the 
envelope are 106 Sn and 104 > 106 ln at the end of the first burst, this switches to 64 Zn, 68 Se, and 
32 S (a similar trend is also reported by Fisker et al. 2008). In this work, the mass-averaged 
composition at the end of the first burst computed for model 3 (see Table 8) is dominated 
(aside from some residual H and 4 He) by the presence of 104 Pd (0.05, by mass) and 105 Ag 
(0.08), while X( 64 Zn) ~ 0.04. But the peak at the end of the abundance distribution (see 
Fig. I3TJJ) increases with subsequent bursts up to a plateau value, which indicates that these 
heavy nuclei are still produced in similar quantities. This is very different to the results 
reported by Woosley et al. (2004). Indeed, at the end of the fifth burst, the abundance 
pattern, shows still a significant presence of heavy species (i.e., X( 105 Ag) ~ 0.1, X( 104 Pd) 
~ 0.08, and X( 94 Tc) ~ 0.05), together with a simultaneous increase in the abundances of 
intermediate-mass elements, such as 60 Ni (0.06), 64 Zn (0.09), 68 Ge (0.07), or 72 Se (0.04). It is 
finally worth noting that, in agreement with all previous hydrodynamic studies, both models 
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1 and 3 yield very small post-burst abundances of 12 C, below the threshold amount required 
to power superbursts. Even though only a few bursts have been computed for these models, 
they already show a trend on the amount of 12 C that may be expected after many more 
bursts. 

Finally, it is also worth mentioning that large differences exist between the hydrody- 
namic simulations reported here (see also Woosley et al. 2004, and Fisker et al. 2008) and 
those based on one-zone models (i.e., Schatz et al. 1999, 2001) as regards the shape of the 
light curve accompanying the bursting episodes (the primary difference being the presence 
of a long-lasting plateau in the latter). 

The origin of the discrepancies reported is not totally clear and would require additional 
hydrodynamic studies. Notice, however, that the local surface gravity of our model is some- 
what smaller than that adopted in the abovementioned works: whereas a 10 km radius is 
assigned to the 1.4 M Q neutron star in Woosley et al. (2004) (g = 1.86 x 10 14 cm s -2 ), the in- 
tegration of the neutron star structure from the core to its surface, in hydrostatic equilibrium, 
yielded 13.1 km (14.3 km, after general relativity corrections are introduced; see Subsection 
5.1), for our 1.4 M Q neutron star (corresponding to a surface gravity of g = 1.08 x 10 14 cm 
s~ 2 ); in turn, the calculations reported by Fisker et al. (2008), in a general relativity frame- 
work, relied on a 11 km (1.4 M ) neutron star, for which g = 1.53 x 10 14 cm s -2 . Although 
XRB properties depend weakly upon the neutron star mass (or surface gravity), part of the 
differences outlined between the three studies can be attributed to the combined effect of the 
adopted neutron star size (surface gravity) and to differences in the input physics (i.e., nu- 
clear reaction network, opacities, treatment of convection). In particular, the use of Iben's 
opacities may have some effect on the peak luminosities achieved since the larger OPAL 
opacities will likely decrease the amount of energy radiated away from the star. Moreover, 
the inclusion of semiconvection and thermohaline mixing would have a minor effect in the 
properties of the explosions, likely affecting the appearance of marginal convective transport 
between bursts (see Woosley et al. 2004, Fisker et al. 2008). It is however worth noting that 
the convective pattern shown in Figs. [18] & [32] is similar to those reported in previous work: 
namely that convection sets in as soon as superadiabatic gradients are established in the 
envelope, following the early stages of the TNR and the corresponding rise in temperature; 
it reaches the surface and begins to recede before the observed burst properly commences, 
shutting off thereafter (Woosley et al. 2004, Fisker et al. 2008). 

The potential impact of XRB nucleosynthesis on Galactic abundances is still a mat- 
ter of debate. Matter accreted onto a neutron star of mass M and radius R releases 
GMrrip/R ~ 200 MeV nucleon -1 , whereas only a few MeV nucleon -1 are released from 
thermonuclear fusion. Thus ejection from a neutron star is unlikely. However, it has been 
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suggested that radiation- driven winds during photospheric radius expansion may lead to ejec- 
tion of a tiny fraction of the envelope (containing nuclear processed material; see Weinberg et 
al. 2006; MacAlpine et al. 2007). Indeed, XRBs have been proposed as a possible source of 
the light p-nuclei 92,94 Mo and 96 ' 98 Ru (Schatz et al. 1998, 2001). No matter is ejected in any 
of the models reported in this work, a result fully independent of the adopted resolution and 
in agreement with all previous hydrodynamic simulations (Woosley et al. 2004, Fisker et al. 
2008). Moreover, it is worth noting that, as shown in Figs. [17]&|30j the abundances of many 
species synthesized during the bursts decrease remarkably towards the outer envelope layers, 
because of inefficient convective transport (see Figs. [18]&[32]). To assess the possible contri- 
bution to the Galactic abundances, one has to rely on the abundances of the outer envelope 
layers (the only ones that have a chance to be ejected by radiation-driven winds). This shows 
the limitations posed by one-zone nucleosynthesis calculations, in which the chemical species 
synthesized in the innermost layers are, by construction, assumed to represent the whole 
(chemically homogeneous) envelope. The mass fractions of these p-nuclei, obtained in model 
3, drop by more than an order of magnitude in the outer envelope layers (as compared with 
the values achieved at the innermost envelope; see Fig. 134")) ; the resulting overproduction 
factors, / ~ 10 6 , are several orders of magnitude smaller than those required to account for 
the origin of these problematic nuclei (see Weinberg et al. 2006, Bazin et al. 2008), in sharp 
contrast with the results obtained on the basis of one-zone calculations (Schatz et al. 1998, 
2001). 
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Table 1. Summary of the models computed in this work. 



Model Mns (M ) Metallicity Envelope shells Bursts computed 

4 
2 
5 



J A mass-accretion rate of 1.75 x 10 9 M yr l has been adopted for all models. 



1 


1.4 


0.02 


60 


2 


1.4 


0.02 


200 


3 


1.4 


1 x 10~ 3 


60 
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Table 2. 



Mean composition of the envelope (5Q > 10 9 ) 

Model 1. 



at the end of each burst, for 



Nucleus Burst 1 



Burst 2 



Burst 3 



Burst 4 



J H 


1.7 


X 


io- 1 


8.8 


X 


io- 2 


6.3 


X 


io- 2 


3.9 


x IO" 2 


4 He 


3.1 


X 


io- 1 


2.0 


X 


io- 1 


1.4 


X 


io- 1 


1.0 


x IO" 1 


12 C 


1.7 


X 


io- 2 


2.5 


X 


io- 2 


2.0 


X 


io- 2 


1.9 


x IO" 2 


13 C 


8.4 


X 


10~ 5 


1.4 


X 


io- 4 


1.3 


X 


io- 4 


1.1 


x IO" 4 


14 N 


2.1 


X 


io- 3 


1.5 


X 


io- 3 


1.1 


X 


io- 3 


7.6 


x IO" 4 


15 N 


2.6 


X 


IO" 3 


2.0 


X 


io- 3 


1.4 


X 


io- 3 


9.6 


x IO" 4 


16Q 


5.7 


X 


io- 4 


5.2 


X 


IO" 4 


3.3 


X 


io- 4 


2.9 


x IO" 4 


17 


3.7 


X 


io- 6 


6.3 


X 


io- 6 


6.8 


X 


io- 6 


1.2 


x IO- 5 


18Q 


6.7 


X 


io- 5 


2.3 


X 


io- 5 


1.7 


X 


io- 5 


6.1 


x IO- 6 


18p 


6.2 


X 


io- 5 


5.4 


X 


io- 5 


3.1 


X 


io- 5 


2.5 


x IO- 5 


19p 


2.1 


X 


io- 4 


1.8 


X 


io- 4 


1.1 


X 


io- 4 


7.3 


x IO- 5 


20 Ne 


6.3 


X 


io- 4 


5.4 


X 


io- 4 


3.5 


X 


io- 4 


4.0 


x IO" 4 


21 Ne 


2.0 


X 


io- 5 


9.5 


X 


io- 6 


6.7 


X 


io- 6 


6.2 


x IO- 6 


22 Ne 


9.9 


X 


io- 5 


4.2 


X 


io- 5 


3.2 


X 


io- 5 


2.9 


x IO- 5 


22 Na 


3.3 


X 


io- 3 


1.9 


X 


io- 3 


1.2 


X 


io- 3 


9.2 


x IO- 4 


23 Na 


3.7 


X 


IO" 4 


2.1 


X 


10~ 4 


1.2 


X 


10~ 4 


1.3 


x IO" 4 


24 Mg 


1.7 


X 


io- 3 


1.4 


X 


io- 3 


9.1 


X 


io- 4 


1.3 


x IO- 3 


25 Mg 


2.6 


X 


io- 3 


1.7 


X 


io- 3 


1.2 


X 


io- 3 


1.5 


x IO- 3 


26 Mg 


1.8 


X 


io- 3 


2.5 


X 


10~ 3 


1.7 


X 


io- 3 


1.1 


x IO- 3 


26 A1S 


2.6 


X 


io- 4 


9.8 


X 


io- 5 


5.7 


X 


io- 5 


1.1 


x IO" 4 


27 Al 


1.8 


X 


io- 3 


1.5 


X 


io- 3 


8.8 


X 


io- 4 


6.7 


x IO" 4 


28 Si 


1.2 


X 


io- 3 


5.3 


X 


io- 3 


7.0 


X 


io- 3 


1.3 


x IO" 2 


29 Si 


3.1 


X 


io- 4 


8.6 


X 


io- 4 


4.6 


X 


io- 4 


6.9 


x IO" 4 


30 Si 


3.5 


X 


io- 3 


5.6 


X 


io- 3 


3.8 


X 


io- 3 


3.0 


x IO- 3 


31p 


5.1 


X 


10~ 4 


1.4 


X 


io- 3 


1.7 


X 


io- 3 


1.8 


x IO- 3 


32 s 


3.8 


X 


10~ 4 


2.9 


X 


io- 2 


5.8 


X 


io- 2 


9.0 


x IO- 2 


33 s 


2.9 


X 


IO" 4 


3.1 


X 


io- 3 


4.3 


X 


io- 3 


5.1 


x IO- 3 


34 s 


3.2 


X 


io- 3 


1.1 


X 


io- 2 


1.6 


X 


io- 2 


1.8 


x IO" 2 


35Q 


1.0 


X 


io- 3 


4.8 


X 


io- 3 


1.0 


X 


io- 2 


1.0 


x IO" 2 


36 Ar 


4.2 


X 


io- 4 


4.0 


X 


io- 3 


9.2 


X 


io- 3 


1.0 


x IO" 2 


37 C1 


3.6 


X 


io- 7 


6.4 


X 


io- 7 


1.8 


X 


io- 6 


4.1 


x IO- 6 


37 Ar 


2.3 


X 


io- 4 


5.7 


X 


io- 4 


1.0 


X 


io- 3 


1.1 


x IO- 3 


38 Ar 


1.9 


X 


io- 3 


4.9 


X 


io- 3 


8.5 


X 


io- 3 


8.9 


x IO- 3 
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Table 2 — Continued 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 



39 K 

40 Ca 



3.8 x 10~ 3 7.8 x lO- 3 1.2 x 10~ 2 1.3 x 10~ 2 

3.2 x KT 3 4.3 x 10~ 3 5.6 x HT 3 5.3 x 10~ 3 
41 K - 9.3 x 10" 9 4.8 x 1(T 8 7.6 x 1(T 8 

Ca 6.5 x 10~ 5 6.9 x KT 5 6.0 x KT 5 7.2 x 10~ 5 

5.5 x 10~ 4 1.2 x 10~ 3 1.7 x 1(T 3 1.9 x 10~ 3 

1.8 x 10~ 4 3.6 x 1(T 4 7.7 x 10~ 4 1.2 x 10~ 3 

4.3 x 1(T 4 9.5 x 10~ 4 1.2 x 10~ 3 9.7 x 10~ 4 

4.3 x 1(T 8 3.7xl0~ 6 1.8 x 10~ 5 3.2 x 1(T 5 
2.5 x 10~ 7 1.9 x 10~ 5 3.6 x 1(T 5 2.9 x 10~ 5 
5.8 x 10~ 4 7.2 x 1(T 4 7.3 x 10~ 4 7.5 x 1(T 4 
8.2 x 10~ 5 7.9 x 1(T 5 9.3 x 1(T 5 1.2 x 10~ 4 
1.5 x 1(T 4 1.7 xlO~ 4 1.1 x 1(T 4 1.0 x 10 -4 

9.5 x 10~ 4 1.9 x 10~ 3 2.4 x 10~ 3 2.7 x lO -3 

6.6 x 10~ 4 8.7 x 10~ 4 9.2 x 10~ 4 9.5 x 10~ 4 

2.1 x 10~ 7 5.6 x 10~ 7 1.6 x 10- 6 3.4 x 10~ 6 

1.2 x 10~ 4 1.4 x 10~ 4 2.4 x KT 4 3.7 x 10~ 4 
1.8 x 10~ 3 2.1 x lO" 3 2.0 x 10~ 3 2.0 x 10~ 3 

1.4 x 10~ 7 4.4 x 10~ 7 6.9 x 10~ 7 8.9 x 10~ 7 

1.5 x 10~ 3 1.5 x 10~ 3 1.2 x 10~ 3 1.4 x 10~ 3 
1.5 x 10~ 3 2.3 x 10~ 3 2.7 x 10~ 3 3.0 x 10~ 3 
3.4 x KT 6 6.4 x 10~ 6 1.6 x 10~ 5 3.0 x 10~ 5 

Cr 3.2 x 10~ 3 4.5 x 10~ 3 5.1 x 10~ 3 5.8 x 10~ 3 



4ir<a p. k xx in-5 

42 Ca 
43 Ca 
43 Sc 
44 Ca 
44 Sc 

44 Ti 

45 Sc 

45 Ti 
46 Ti 
47 Ti 
48 Ti 

48y 1 o sx 1 n- z; i a nx 1 n- z! o a nx 1 n- z! q 7 nx 1 n- zi 

48 Cr 
49 T j 
49y 
50 Cr 

5i y q a nx m-<> a a nx in-i' 1 a nx in-5 q n nx in-5 

51 

52 Cr 

52 Mn 

52p e 

53 Mn 

54 Fe 

55 Mn 
55 Fe 

55 Co 
56 Fe 
56 Co 

56 Ni 



1.5 x 10~ 5 2.2 x 10~ 5 4.9 x 10~ 5 9.4 x 10~ 5 
3.4 x lO" 3 2.3 x 10~ 3 2.7 x 10~ 3 4.0 x 10~ 3 
2.0 x lO" 2 1.5 x 10~ 2 9.9 x 10~ 3 1.2 x 10~ 2 

1.2 x lO" 3 1.1 x 10~ 3 8.6 x 10~ 4 1.1 x 10~ 3 
1.0 x 10~ 3 1.4 x 10~ 3 1.6 x 10~ 3 1.8 x 10~ 3 

6.6 x 10~ 9 2.3 x 10~ 8 7.1 x 10~ 8 1.6 x 10~ 7 

2.3 x 10~ 4 3.1 x 10~ 4 6.4 x KT 4 1.0 x 10~ 3 
2.9 x 10~ 3 3.8 x 10~ 3 4.0 x KT 3 4.3 x 10~ 3 
8.0 x 10~ 8 5.4xl0~ 7 1.7xl0~ 6 3.5 x 10~ 6 
2.2 x lO" 4 5.3 x 10~ 4 1.1 x 10~ 3 1.8 x 10~ 3 

2.4 x 10- 2 4.5 x 10~ 2 5.2 x 10~ 2 5.9 x 10~ 2 
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Table 2 — Continued 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 



2.6 x 1(T 8 1.6 x 10~ 7 2.0 x 10~ 7 2.5 x 10~ 7 
2.5 x 10~ 4 1.2 x 10" 4 1.4 x 10~ 4 2.0 x 10~ 4 
6.5 x 10~ 3 3.9 x 10~ 3 2.7 x 10~ 3 2.9 x 10~ 3 

4.7 x 10~ 3 2.7 x 10~ 3 1.9 x 10" 3 2.2 x 1(T 3 

7.3 x 10~ 3 4.3 x 10~ 3 3.3 x 10~ 3 3.8 x 10" 3 
3.2 x lO" 1 3.3 x 10" 1 3.1 x 10" 1 3.1 x 10" 1 
3.9 x 10" 3 1.7 x 10" 3 1.9 x 10" 3 2.4 x 10" 3 
8.2 x lO" 3 6.0 x 10~ 3 3.5 x 10" 3 3.8 x 10" 3 

5.0 x 10~ 4 3.1 x 10" 4 4.4 x 10' 4 5.8 x 10~ 4 
3.2 x 10~ 3 2.5 x 10" 3 1.4 x 10" 3 1.5 x 10~ 3 

4.1 x 10~ 3 4.4 x 10~ 3 3.1 x 10" 3 3.3 x 10" 3 

3.4 x 10~ 2 1.0 x lO" 1 1.3 x 10" 1 1.3 x 10" 1 

2.5 x 10~ 7 2.9 x 10~ 7 5.1 x 10~ 7 9.2 x 10~ 7 
1.4 x 10~ 3 2.5 x 10~ 3 1.9 x 10~ 3 2.2 x 10~ 3 

3.2 x 10~ 5 5.0 x 10- 5 1.3 x 10~ 4 2.4 x 10~ 4 
4.0 x 10^ 4 5.1 x 10~ 4 5.5 x 10^ 4 6.3 x 10~ 4 

5.4 x 10~ 4 1.3 x 10~ 3 7.5 x 10~ 4 8.7 x 10~ 4 
4.7 x 10~ 4 1.0 x 10~ 3 9.7 x 10~ 4 1.1 x 10~ 3 
2.9 x 10~ 3 2.6 x 10~ 2 4.4 x 10~ 2 4.5 x 10^ 2 

2.2 x 10~ 4 1.3 x 10~ 3 1.5 x 10~ 3 1.8 x 10~ 3 

7.3 x 10~ 5 4.1 x 10~ 4 4.6 x HT 4 6.0 x 10" 4 
7.7 x 10~ 5 4.7 x 10~ 4 6.0 x 10~ 4 7.7 x 10~ 4 

2.5 x 10~ 4 5.7 x 10~ 3 1.4 x 10~ 2 1.6 x 10" 2 
3.0 x lO" 5 4.2 x 10~ 4 7.6 x 10~ 4 1.0 x 10~ 3 

1.3 x 10~ 5 1.8 x 10~ 4 3.3 x 10~ 4 4.8 x 10" 4 
1.2 x lO" 5 1.8 x 10~ 4 3.7xl0~ 4 5.3 x 10~ 4 

Kr 3.0 x 10~ 5 1.4 x 10~ 3 4.5 x 10~ 3 5.6 x 10~ 3 

4.5 x 10~ 6 1.5 x 10~ 4 4.1 x 10~ 4 6.1 x 10~ 4 
2.7 x 10~ 6 9.2 x 10~ 5 2.5 x 10~ 4 4.1 x 10~ 4 

1.6 x 10~ 6 5.9 x 10~ 5 1.7 xl(T 4 2.7 x 10~ 4 
3.9 x 10~ 6 4.0 x 10~ 4 1.6 x 10~ 3 2.2 x 10~ 3 

8.4 x 10~ 7 6.8 x 10~ 5 2.7 x 10~ 4 4.3 x 10~ 4 
7.9 x 10~ 7 8.2 x lO" 5 3.2 x 10~ 4 5.6 x 10~ 4 



57 Fe 

57 Co 

57 Ni 
58 Ni 

60 Ni 
61 Ni 

61 Cu 

62 M 

62 Zn 
63 Cu 
64 Zn 
65 Cu 
65 Zn 
66 Zn 
66 Ga 
66 Ge 
67 Ga 
68 Ge 
69 Ge 
70 Ge 
71 As 
72 Se 
73 Se 
74 Se 

75 Br 



76 

77 Kr 
78 Kr 
79 Kr 

80 Sr 
81 Rb 

82g r 
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Table 2 — Continued 



Nucleus Burst 1 



Burst 2 



Burst 3 



Burst 4 



83 Sr 
84 Sr 

85y 

86 Zr 
87 Zr 

88 Zr 

89 Nb 

90 Mo 
91 Nb 
92 Mo 

93 Tc 
94 Tc 

95 Ru 

96 Ru 

97 Ru 

98 Ru 

"Rh 
ioo pd 



3.0 x 10- 7 4.3 x 10- 5 
1.9 x 10- 7 3.8 x 10- 5 

1.1 x 10- 7 4.0 x 10- 5 
7.1 x 10- 8 4.4 xlO- 5 
2.7 x 10- 8 3.6 x 10- 5 
5.6 x 10- 9 1.8 x 10- 5 
4.0 x 10- 9 3.1 x 10- 5 

1.1 x 10- 5 
5.4 x 10- 6 
2.1 x 10- 6 

1.3 x 10- 6 
6.9 x 10- 7 
1.1 x 10- 7 

9.4 x 10- 9 
3.9 x 10- 9 



1.9 x 


io- 4 


3.2 


x IO" 4 


1.4 x 


io- 4 


2.6 


x IO" 4 


2.1 x 


io- 4 


3.5 


x IO" 4 


2.3 x 


io- 4 


4.6 


x IO" 4 


1.8 x 


io- 4 


3.9 


x IO- 4 


1.4 x 


io- 4 


2.6 


x IO" 4 


2.9 x 


io- 4 


5.8 


x IO" 4 


1.1 x 


io- 4 


2.2 


x IO" 4 


5.3 x 


io- 5 


1.1 


x IO" 4 


1.9 x 


io- 5 


4.5 


x IO- 5 


1.3 x 


io- 5 


4.5 


x IO- 5 


7.6 x 


io- 6 


3.8 


x IO- 5 


1.1 x 


io- 6 


4.7 


x IO- 6 


1.1 x 


io- 7 


7.6 


x IO" 7 


5.1 x 


io- 8 


3.9 


x IO" 7 


1.3 x 


io- 8 


9.2 


x IO- 8 


2.4 x 


io- 9 


1.5 


x IO- 8 


. 




2.2 


x IO" 9 



Table 3. Summary of burst properties for Model 1. 



Burst T peak (K) t(T P eak) (s) r rec (hr) L peak (L ) r .oi (s) a 



1 


1.06 x IO 9 


21192 


5.9 


9.7 x IO 4 


75.8 


60 


2 


1.15 x IO 9 


44342 


6.4 


1.7 x IO 5 


62.3 


40 


3 


1.26 x IO 9 


62137 


4.9 


2.1 x IO 5 


55.4 


34 


4 


1.12 x IO 9 


80568 


5.1 


1.2 x IO 5 


75.7 


36 
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Table 4. Summary of burst properties for Model 2. 



Burst 


J-peafc \^-) 


tijpeak) (s) 


Tree (hr) 


Lpeafc (L©) 


TO.01 (s) 


ct 


1 

2 


1.05 x 10 9 
1.20 x 10 9 


21189 
37783 


5.9 
4.6 


9.0 x 10 4 
1.5 x 10 5 


59.2 
73.9 


62 
31 
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Table 5. 



Mean composition of the envelope (5Q > 10 9 ) 

Model 2. 



at the end of each burst, for 



Nucleus Burst 1 



Burst 2 



J H 


2.3 


X 


io- 1 


1.0 


x IO" 1 


4 He 


3.0 


X 


io- 1 


1.9 


x IO" 1 


12 C 


1.8 


X 


io- 2 


2.2 


x IO" 2 


13 C 


6.8 


X 


10~ 5 


6.2 


x IO" 5 


14 N 


2.0 


X 


io- 3 


1.4 


x IO- 3 


15 N 


2.6 


X 


io- 3 


2.0 


x IO- 3 


16Q 


5.3 


X 


io- 4 


3.9 


x IO" 4 


17 


4.8 


X 


IO" 6 


3.0 


x IO" 6 


18Q 


6.3 


X 


io- 5 


2.8 


x IO" 5 


18p 


2.1 


X 


io- 4 


5.7 


x IO" 5 


19p 


3.6 


X 


io- 4 


1.7 


x IO" 4 


20 Ne 


5.4 


X 


io- 4 


4.1 


x IO" 4 


21 Ne 


2.2 


X 


io- 5 


7.3 


x IO" 6 


22 Ne 


7.7 


X 


io- 5 


4.4 


x IO" 5 


22 Na 


3.6 


X 


io- 3 


1.6 


x IO- 3 


23 Na 


2.3 


X 


io- 4 


1.7 


x IO" 4 


24 Mg 


1.2 


X 


io- 3 


1.3 


x IO- 3 


25 Mg 


1.8 


X 


io- 3 


1.4 


x IO- 3 


26 Mg 


3.1 


X 


io- 3 


1.6 


x IO- 3 


26 A1S 


1.0 


X 


io- 4 


1.1 


x IO" 4 


27 Al 


1.6 


X 


io- 3 


1.1 


x IO- 3 


28 Si 


1.2 


X 


io- 3 


7.5 


x IO- 3 


29 Si 


4.6 


X 


io- 4 


6.7 


x IO" 4 


30 Si 


4.7 


X 


10~ 3 


4.6 


x IO- 3 


31p 


6.1 


X 


10~ 4 


1.4 


x IO- 3 


32 s 


4.7 


X 


io- 4 


2.7 


x IO" 2 


33 s 


4.0 


X 


IO" 4 


2.4 


x IO- 3 


34 s 


4.3 


X 


io- 3 


9.6 


x IO- 3 


35Q 


1.1 


X 


io- 3 


3.8 


x IO- 3 


36 Ar 


5.8 


X 


io- 4 


2.9 


x IO- 3 


37 C1 


1.6 


X 


io- 7 


5.3 


x IO" 7 


37 Ar 


2.9 


X 


io- 4 


4.0 


x IO" 4 


38 Ar 


2.6 


X 


io- 3 


5.2 


x IO- 3 
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Table 5 — Continued 



Nucleus Burst 1 Burst 2 



4.9 x 10~ 3 8.7 x 10~ 3 

3.1 x 10~ 3 3.4xl0~ 3 

41 K - 8.7 x KT 9 

41 Ca 6.8 x KT 5 4.5 x 10~ 5 

42 Ca 7.3 x 10~ 4 1.2 x lO" 3 

43 Ca 9.6 x 1(T 5 4.3 x 10~ 4 

43 Sc 7.2 x 1(T 4 1.2 x 1(T 3 

44 Ca 5.8 x 10~ 9 2.3 x 1(T 6 

9.6 x 10- 8 1.1 x 10~ 5 

5.6 x 1(T 4 4.8 x 10~ 4 

3.5 x 1(T 5 4.9 x 1(T 5 

2.1 x 10~ 4 1.2 x 1(T 4 

1.2 x 10~ 3 1.9 x 10~ 3 
6.8 x 10~ 4 6.3 x 10~ 4 

2.8 x 1(T 8 2.5 x 10~ 7 
V 4.4 x 1(T 5 9.4 x 1(T 5 

1.9 x 10" 3 1.7 x 1CT 3 
2.4 x 10~ 8 2.9 x 10~ 7 
1.4 x 10~ 3 1.2 x 10~ 3 

1.6 x 1(T 3 2.2 x 10~ 3 
6.4 x 10" 7 5.5 x 1(T 6 
3.4 x 10- 3 4.9 x 10- 3 
1.6 x 10~ 6 8.1 x 10~ 6 
1.0 x 10" 3 1.5 x 10- 3 

Fe 1.7 x 10~ 2 1.3 x 1(T 2 

Mn 9.9 x 1(T 4 8.4 x 10~ 4 

Fe 9.2 x 10~ 4 1.2 x 10" 3 



39 K 

40 Ca 



44 Sc 

44 Ti 

45 Sc 

45 Ti 
46 Ti 
47 Ti 
48 Ti 

18 
48 Cr 
49 T j 

49y 
50 Cr 

51y 

51 Cr 

52 Cr 

52 Mn 

52 
53 

51 



55 Mn - 1.3 x 1CT 8 

55 Fe 

55 Co 2.8 x 10" 3 3.7 x 10" 3 

56 Fe 

56 Co 

56 Ni 



7.9 x 10~ 5 2.6 x 10~ 4 

2.8 x 10" 3 3.7 x 10" 3 

1.4 x 10~ 8 3.6 x 1(T 7 

1.1 x lO" 4 5.0 x 10~ 4 

3.2 x 10~ 2 5.0 x 10~ 2 
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Table 5 — Continued 



Nucleus Burst 1 Burst 2 



2.6 x 10~ 9 6.6 x 10- 8 

6.9 x 10~ 5 7.4 x 10~ 5 

5.1 x 10~ 3 3.4xl0- 3 

3.4 x 1(T 3 2.3 x 10~ 3 

4.8 x 1(T 3 3.7 x 1(T 3 

2.9 x 10" 1 3.2 x KT 1 
9.1 x 10~ 4 1.4xlO~ 3 
6.3 x 1(T 3 5.8 x lO" 3 
9.9 x 10~ 5 1.9 x 10" 4 

Zn 1.8 x 10~ 3 2.1 x 1(T 3 

Cu 2.2 x 1(T 3 3.2 x 10~ 3 

Zn 2.3 x 1(T 2 1.2 x 10" 1 

Cu 2.4 x 1(T 8 2.7 x 10~ 7 

Zn 5.4 x 10~ 4 3.1 x 10~ 3 

Zn 1.7 x 10- 6 2.9 x 1(T 5 

Ga 6.4 x 10~ 5 6.2 x 10~ 4 

Ge 2.8 x KT 4 1.7 x 10~ 3 

Ga 1.8 x 1(T 4 1.2 x 10~ 3 

Ge 1.6 x 10~ 3 2.8 x 10~ 2 

Ge 7.3 x 10~ 5 1.8 x 10" 3 

Ge 2.0 x 10~ 5 5.9 x 10~ 4 

As 2.4 x 10~ 5 6.2 x 10~ 4 

; Se 1.2 x KT 4 5.8 x lO -3 

; Se 8.0 x 10~ 6 6.1 x 10~ 4 

Se 3.2 x 10~ 6 2.7 x 10~ 4 

Br 3.3 x 10~ 6 2.7 x 10~ 4 

Kr 1.2 x 10~ 5 1.6 x 10" 3 

Kr 1.2 x 10~ 6 2.3 x 10~ 4 

Kr 6.7 x 10~ 7 1.5 x 10~ 4 

Kr 4.0 x 10~ 7 9.2 x 10~ 5 

'Sr 1.5 x 10~ 6 4.9 x 10~ 4 

Rb 2.3 x 10~ 7 1.1 x 10~ 4 

! Sr 2.2 x 10~ 7 1.4 x 10~ 4 
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Table 5 — Continued 



Nucleus Burst 1 Burst 2 



83 Sr 9.4 x 10- 8 7.0 x lO" 5 

84 Sr 6.4 x 1(T 8 6.3 x 10~ 5 

85 Y 4.5 x 1(T 8 6.4 x 1(T 5 

86 Zr 3.1 x 10~ 8 7.8 x 1(T 5 

87 Zr 1.4 x 10- 8 6.9 x 10~ 5 

88 Zr 3.5 x 10~ 9 3.3 x 1(T 5 

89 Nb 2.5 x 10~ 9 5.7 x 10~ 5 

90 Mo 2.4 x 1(T 5 

91 Nb - 1.3 x 10- 5 

92 Mo - 7.3 x 1(T 6 

93 Tc 6.9 x KT 6 

94 Tc 5.0 x 10~ 6 

95 Ru 9.2 x 10~ 7 

96 Ru 9.1 x 10~ 8 

97 Ru 3.7 x 10- 8 

98 Ru 9.1 x 10~ 9 

"Rh - 1.8 x 10~ 9 
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Table 6. Properties of the last burst computed in Models 1 & 3. 



Model 1 



Model 3 



Pmax,base ^g,-CIii ) 

Pmax,base (dyn.Cm~ 2 ) 

Pmax,ign ^g.CUL J 

Pmax,ign (dyn.cm~ 2 ) 

Tacc (hr) 

Am acc (M ) 

1- peak \^-) 
'-'peak (h©) 
t-^Zmax \"*) 

a 



1.3 x 10 6 2.6 x 10 6 

4.2 x 10 22 1.1 x 10 23 

5.4 x 10 5 1.0 x 10 6 

1.3 x 10 22 2.9 x 10 22 
5.1 8.8 

1.0 x 10~ 12 1.8 x 10~ 12 

1.1 x 10 9 1.3 x 10 9 

1.2 x 10 5 1.0 x 10 5 
40 44 

36 30 



a Pmax,base and P ma x,base are the maximum density and pressure achieved at the base of the 
envelope, whereas p m ax,i gn and P ma x,ign correspond to the maximum values attained at the 
ignition shell (defined as the first shell that reaches T > 4 x 10 8 K, for this burst). 



Table 7. Summary of burst properties for Model 3. 



Burst T peak (K) t(T peafc ) (s) r rec (hr) L peak (L ) r .oi (s) a 



1 


1.40 x 10 9 


65110 


18.1 


1.0 x 10 5 


423 


34 


2 


1.39 x 10 9 


98879 


9.4 


1.1 x 10 5 


296 


24 


3 


1.32 x 10 9 


130816 


8.9 


9.8 x 10 4 


281 


24 


4 


1.30 x 10 9 


162777 


8.9 


1.0 x 10 5 


252 


27 


5 


1.26 x 10 9 


194266 


8.8 


1.0 x 10 5 


250 


30 
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Table 8. Mean composition of the envelope (Xj > 10 9 ) 

Model 3. 



at the end of each burst, for 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 Burst 5 



X H 


1.8 


X 


io- 1 


8.6 


X 


io- 2 


5.5 


X 


io- 2 


4.5 


X 


io- 2 


3.7 


x IO" 2 


4 He 


8.4 


X 


io- 2 


5.8 


X 


io- 2 


4.4 


X 


io- 2 


3.8 


X 


io- 2 


3.3 


x IO" 2 


12 C 


7.7 


X 


io- 4 


1.7 


X 


io- 3 


1.8 


X 


io-- 3 


1.7 


X 


io-- 3 


1.5 


x IO- 3 


13 C 


4.2 


X 


io- 6 


8.9 


X 


io- 6 


1.4 


X 


io- 5 


2.3 


X 


io- 5 


2.7 


x IO- 5 


14 N 


2.8 


X 


10" 4 


2.1 


X 


io- 4 


1.7 


X 


10~ 4 


1.5 


X 


10~ 4 


1.4 


x IO" 4 


15 N 


4.9 


X 


io- 4 


3.6 


X 


10~ 4 


2.9 


X 


io- 4 


2.4 


X 


IO" 4 


2.2 


x IO" 4 


16Q 


1.4 


X 


io- 5 


2.9 


X 


io- 5 


4.1 


X 


io- 5 


3.7 


X 


io- 5 


3.8 


x IO- 5 


17 


4.8 


X 


io- 7 


6.8 


X 


io- 7 


2.1 


X 


io- 6 


1.3 


X 


io- 6 


2.0 


x IO- 6 


18Q 








3.7 


X 


io- 9 


9.4 


X 


io- 9 


9.5 


X 


io- 9 


8.1 


x IO" 9 


18p 


6.4 


X 


io- 9 


1.3 


X 


io- 8 


2.8 


X 


io- 8 


2.4 


X 


io- 8 


2.3 


x IO- 8 


19p 


8.7 


X 


io- 8 


1.6 


X 


io- 7 


1.7 


X 


io- 7 


2.0 


X 


io- 7 


8.0 


x IO- 8 


20 Ne 


1.5 


X 


io- 5 


2.9 


X 


io- 5 


3.5 


X 


io- 5 


3.0 


X 


io- 5 


3.1 


x IO- 5 


21 Ne 


1.9 


X 


io- 8 


4.4 


X 


io- 8 


4.9 


X 


io- 8 


5.2 


X 


io- 8 


4.4 


x IO- 8 


22 Ne 


2.9 


X 


io- 7 


6.6 


X 


io- 7 


4.3 


X 


io- 7 


1.1 


X 


io- 6 


7.5 


x IO" 7 


22 Na 


1.1 


X 


io- 5 


2.4 


X 


io- 5 


2.0 


X 


io- 5 


2.3 


X 


io- 5 


1.9 


x IO- 5 


23 Na 


2.6 


X 


io- 6 


6.5 


X 


io- 6 


6.5 


X 


io- 6 


6.5 


X 


io- 6 


5.7 


x IO- 6 


24 Mg 


8.2 


X 


io- 5 


1.4 


X 


io- 4 


1.4 


X 


io- 4 


1.2 


X 


10~ 4 


1.1 


x IO" 4 


25 Mg 


4.9 


X 


io- 5 


7.3 


X 


io- 5 


7.4 


X 


io- 5 


7.2 


X 


io- 5 


6.6 


x IO- 5 


26 Mg 


1.1 


X 


io- 5 


1.9 


X 


io- 5 


1.8 


X 


io- 5 


2.2 


X 


io- 5 


2.0 


x IO- 5 


26 A1S 


1.2 


X 


io- 5 


1.4 


X 


io- 5 


1.3 


X 


io- 5 


1.2 


X 


io- 5 


1.1 


x IO- 5 


27 Al 


4.1 


X 


io- 5 


6.3 


X 


io- 5 


6.5 


X 


io- 5 


6.8 


X 


io- 5 


6.2 


x IO- 5 


28 Si 


2.6 


X 


io- 4 


6.2 


X 


io- 4 


5.2 


X 


io- 4 


3.9 


X 


io- 4 


3.6 


x IO" 4 


29 Si 


8.6 


X 


io- 6 


1.1 


X 


io- 5 


9.2 


X 


io- 6 


1.3 


X 


io- 5 


1.1 


x IO- 5 


30 Si 


5.5 


X 


io- 5 


8.8 


X 


io- 5 


8.0 


X 


io- 5 


1.1 


X 


io- 4 


9.9 


x IO- 5 


31p 


2.3 


X 


io- 5 


3.6 


X 


io- 5 


2.6 


X 


io- 5 


3.2 


X 


io- 5 


2.6 


x IO- 5 


32 s 


3.7 


X 


io- 5 


2.1 


X 


io- 3 


4.9 


X 


io- 3 


5.8 


X 


IO" 3 


7.7 


x IO- 3 


33 s 


5.7 


X 


io- 6 


4.2 


X 


io- 5 


6.1 


X 


io- 5 


5.6 


X 


io- 5 


9.2 


x IO- 5 


34 s 


2.6 


X 


io- 5 


6.7 


X 


io- 5 


7.2 


X 


io- 5 


7.3 


X 


io- 5 


7.5 


x IO- 5 


35Q 


2.9 


X 


io- 5 


1.5 


X 


io- 4 


4.2 


X 


io- 4 


5.0 


X 


io- 4 


5.9 


x IO" 4 


36 Ar 


8.2 


X 


io- 6 


2.9 


X 


io-- 3 


6.8 


X 


io-- 3 


9.9 


X 


io- 3 


1.2 


x IO" 2 


37 C1 


3.5 


X 


io- 9 


2.8 


X 


io- 8 


5.4 


X 


io- 8 


1.5 


X 


io- 7 


2.2 


x IO" 7 


37 Ar 


3.7 


X 


io- 6 


2.6 


X 


io- 5 


5.1 


X 


io- 5 


1.1 


X 


io- 4 


1.6 


x IO" 4 


38 Ar 


1.5 


X 


1Q- 5 


6.4 


X 


1Q- 5 


1.1 


X 


io- 4 


1.7 


X 


io- 4 


2.0 


x IO" 4 
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Table 8 — Continued 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 Burst 5 



39 K 4.7 x 10~ 5 9.4 x 1(T 4 1.4 x 1(T 3 2.1 x 10~ 3 2.5 x 10~ 3 

40 Ca 1.3 x 10~ 4 1.3 x 10~ 3 1.8 x 10~ 3 2.3 x 10~ 3 2.6 x 1(T 3 

41 K - - - - 1.1 x 1(T 9 

41 Ca 2.5 x 10~ 6 6.4 x 1(T 6 7.7 x 10" 6 1.4 x 10~ 5 1.7 x 1(T 5 

42 Ca 3.7 x 10- 6 1.6 x 1(T 5 1.8 x 1(T 5 3.3 x lO" 5 3.5 x 1(T 5 

43 Ca 3.3 x 1(T 6 1.1 x 10~ 5 1.0 x 1(T 5 2.0 x 10~ 5 2.6 x 10~ 5 

1.4 x 10~ 5 4.0 x 10~ 5 3.7 x 10~ 5 5.3 x 10~ 5 5.2 x 10" 



43 Sc 1.4 x 10~ 5 4.0 x 10~ 5 3.7 x KT 5 5.3 x 10~ 5 5.2 x 10~ 5 

44 Ca 4.2 x 10~ 9 5.7 x 10~ 7 1.6 x 10~ 6 5.1 x 10~ 6 

2.8 x 10~ 8 2.2 x 10" 6 3.3 x 10~ 6 9.2 x 10~ 6 

/in., m-5 CO., ia-5 on., 1 n-5 a r, ., i n-5 



Sc 6.0 x 10- 9 



! ~" U.U A1U " Z.O A 1U " Z.Z A 1U ' O.O A ±U " y.z A 1U " 

M Ti 1.5 x 10~ 5 4.9 x 10~ 5 6.3 x 10~ 5 8.0 x 10~ 5 9.0 x 10~ 5 

1.2 x 10~ 6 1.9 x 10~ 6 6.0 x 10~ 6 1.3 x 10~ 5 2.3 x 10~ 5 

An., in-6 en., 1 n-6 10., 1 n-5 1 n ., 1 n-5 00., 1 n-5 



L)U l.Z, A. 1U -L.y A 1U U.U A _LU ±.<J A 1U Zj.O A J.U 

45 Ti 4.0 x 10~ 6 6.0 x 10~ 6 1.8 x 10~ 5 1.9 x 10~ 5 2.3 x 10" 

46 Ti 7.0 x 10~ 6 1.4 x 10~ 5 1.8 x 10~ 5 3.9 x 10~ 5 5.4 x 10' 

47 Ti 1.4 x 10~ 5 3.0 x 10~ 5 3.9 x 10~ 5 5.7 x 10~ 5 6.7 x 10~ 5 

48 Ti 2.6 x 10~ 9 5.8 x 10" 9 4.9 x 10~ 8 2.0 x 10~ 7 5.4 x 10~ 7 

48 V 2.4 x 10~ 6 4.8 x 10~ 6 1.3 x 10~ 5 2.7 x 10~ 5 4.7 x 10~ 5 

48 Cr 5.2 x 10~ 5 1.0 x 10" 4 1.2 x 10~ 4 1.3 x 10' 4 1.3 x 10~ 4 

49 Ti 1.9 x 10~ 9 2.2 x 10~ 8 4.8 x 10~ 8 8.3 x 10~ 8 1.2 x 10~ 7 

49 V 4.7 x 10~ 5 6.3 x 10~ 5 7.9 x lO" 5 8.8 x 10~ 5 9.5 x 10~ 5 

50 Cr 3.9 x 10~ 5 5.7 x 10~ 5 6.9 x 10~ 5 1.0 x 10~ 4 1.3 x 10~ 4 

5.3 x lO" 8 1.1 x 10~ 7 3.2 x 10~ 7 7.9 x 10~ 7 1.7 x 10~ 6 



v o.o x ru i.i x iu o.z x iu i .y x iu i.i x iu 

Cr 1.2 x KT 4 1.9 x 10" 4 2.2 x 10" 4 2.7 x 10" 4 3.1 x 10~ 4 

Cr 5.6 x 10~ 7 7.4 x 10~ 7 4.0 x KT 6 1.3 x 10~ 5 2.8 x 10~ 5 

Mn 2.1 x 10~ 4 2.5 x 10~ 4 4.0 x 10~ 4 6.0 x 10~ 4 7.7 x 10~ 4 

Fe 2.1 x KT 3 2.2 x 10~ 3 2.1 x 10~ 3 1.8 x 10~ 3 1.6 x 10~ 3 

Mn 8.8 x 10~ 5 2.1 x 10~ 4 3.6 x 10~ 4 4.0 x 10~ 4 4.2 x 10~ 4 

Fe 7.7 x 10~ 5 1.2 x 10~ 4 1.5 x 10~ 4 1.8 x 10~ 4 1.9 x 10~ 4 

Mn - - 5.1 x 10~ 9 1.6 x 10~ 8 3.4 x 10~ 8 

'Fe 1.0 x 10~ 5 1.7xl0~ 5 4.5 x 10~ 5 7.7 x 10~ 5 1.1 x 10~ 4 

Co 2.2 x 10~ 4 3.1 x 10~ 4 3.0 x KT 4 3.0 x 10~ 4 2.9 x 10~ 4 

Fe 1.0 x 10~ 9 1.3 x 10~ 8 1.4 x 10~ 7 4.7 x 10~ 7 9.9 x 10~ 7 

Co 4.6 x 10~ 6 2.2 x 10~ 5 8.4 x 10~ 5 1.6 x 10" 4 2.5 x 10" 4 

; Ni 7.8 x 10" 4 2.1 x 10~ 3 2.7 x 10~ 3 3.3 x 10~ 3 3.6 x 10" 3 



66 Ge 
67 Ga 
68 Ge 
69 Ge 
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Table 8 — Continued 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 Burst 5 

57 Fe - 3.3 x 10~ 8 6.2 x lO" 8 1.0 x 10~ 7 1.5 x 10~ 7 

57 Co 



O.O A 1U U.Z A 1U l.U A 1U l.O A 1U 

1.3 x 10~ 5 1.5 x 10~ 5 2.9 x lO" 5 4.8 x 1(T 5 6.6 x 1(T 5 
5.6 x 10~ 4 5.8 x 1(T 4 5.7 x 10" 4 5.1 x 10" 4 4.8 x 1(T 4 
4.1 x 10" 4 4.0 x 10~ 4 4.0 xlO -4 3.7 x 10~ 4 3.6 x 10~ 4 
6.9 x 10~ 4 6.4 x 10~ 4 6.3 x 10~ 4 5.8 x 10~ 4 5.7 x 10~ 4 

2.4 x 10~ 2 3.7 x 10~ 2 4.6 x 10~ 2 5.2 x 10~ 2 5.8 x 10" 2 
-m 1.5 x 10~ 3 2.1 x 10- 3 3.8 x 10~ 3 5.5 x 10~ 3 6.3 x 10~ 3 
L Cu 5.7 x 10~ 3 6.9 x 10~ 3 6.5 x 10~ 3 4.7 x 10~ 3 3.7 x 10~ 3 
2 Ni 3.3 x 10~ 4 5.8 x 10^ 4 1.4 x 10~ 3 1.8 x 10' 3 2.0 x 10~ 3 
2 Zn 3.6 x 10~ 3 4.4 x 10~ 3 3.3 x 10~ 3 2.3 x 10~ 3 1.9 x 10~ 3 
! Cu 4.1 x 10~ 3 3.3 x 10~ 3 2.9 x 10~ 3 2.5 x 10~ 3 2.5 x 10~ 3 
*Zn 4.2 x 10~ 2 5.6 x 10~ 2 7.0 x 10~ 2 8.1 x 10~ 2 9.1 x 10~ 2 
5 Cu 8.1 x 10~ 7 1.2 x 10~ 6 4.9 x 10~ 6 1.1 x 10~ 5 1.7 x 10~ 5 
5 Zn 8.3 x 10~ 3 9.7 x 10~ 3 1.2 x 10~ 2 1.2 x 10~ 2 1.2 x 10~ 2 
3 Zn 1.3 x 10~ 4 2.3 x 10" 4 9.9 x 10~ 4 1.8 x 10~ 3 2.3 x 10~ 3 
5 Ga 2.8 x 10~ 3 2.7 x 10~ 3 3.1 x 10~ 3 2.7 x 10" 3 2.2 x 10~ 3 

6.8 x 10~ 3 4.9 x 10" 3 3.8 x 10~ 3 2.3 x 10~ 3 1.9 x 10~ 3 

4.1 x HT 3 3.4 x 10~ 3 3.6 x 10" 3 3.2 x 10~ 3 3.0 x 10~ 3 

2.8 x 10~ 2 4.5 x 10~ 2 5.6 x 10~ 2 6.5 x 10^ 2 7.1 x 10~ 2 

1.0 x 10~ 2 1.0 x 10~ 3 1.2 x 10~ 2 1.1 x 10~ 2 1.1 x 10~ 2 

70 Ge 6.5 x 10~ 3 5.2 x 10~ 3 5.3 x 10" 3 4.6 x 10~ 3 4.2 x 10~ 3 

71 As 6.1 x 10~ 3 4.8 x 10~ 3 5.0 x 10~ 3 4.4 x 10~ 3 4.2 x 10~ 3 

72 Se 1.7 x 10~ 2 2.8 x 10~ 2 3.4 x KT 2 3.8 x 10~ 2 4.1 x 10~ 2 

73 Se 9.1 x 10~ 3 6.5 x 10~ 3 7.0 x 10~ 3 6.4 x 10~ 3 6.3 x 10~ 3 

74 Se 7.2 x 10~ 3 6.5 x 10~ 3 7.7 x 10~ 3 7.4 x 10~ 3 7.2 x 10~ 3 

75 Br 8.4 x 10~ 3 4.5 x 10~ 3 3.3 x 10~ 3 2.2 x 10~ 3 2.0 x 10~ 3 

76 Kr 1.3 x 10~ 2 1.8 x 10~ 2 2.0 x KT 2 2.2 x 10~ 2 2.3 x 10~ 2 

77 Kr 7.2 x 10~ 3 6.0 x 10~ 3 6.2 x 10~ 3 5.9 x 10~ 3 5.7 x 10~ 3 

78 Kr 8.7 x 10~ 3 7.3 x 10~ 3 7.5 x 10~ 3 6.9 x 10~ 3 6.5 x 10~ 3 

79 Kr 7.0 x 10~ 3 3.5 x 10~ 3 3.0 x 10~ 3 2.5 x 10~ 3 2.4 x 10~ 3 

80 Sr 1.3 x 10~ 2 1.3 x 10~ 2 1.4 x 10~ 2 1.5 x 10~ 2 1.5 x 10~ 2 

81 Rb 8.0 x 10~ 3 6.4 x 10~ 3 6.4 x KT 3 6.0 x 10~ 3 5.8 x 10~ 3 

82 Sr 1.7 x 10~ 2 1.3 x 10~ 2 1.2 x 10~ 2 1.2 x 10~ 2 1.1 x 10~ 2 
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Table 8 — Continued 



Nucleus Burst 1 Burst 2 Burst 3 Burst 4 Burst 5 



83 Sr 1.3 x KT 2 1.2 x 10~ 2 1.4 x 1(T 2 1.3 x 1(T 2 1.3 x 1(T 2 

84 Sr 1.3 x 1(T 2 6.1 x 10~ 3 3.8 x KT 3 2.6 x 10" 3 2.2 x 10" 3 

85 Y 9.8 x 10" 3 7.9 x 10- 3 7.7 x 10" 3 7.2 x 10" 3 6.9 x 10" 3 

86 Zr 1.7 x 1(T 2 1.9 x 10~ 2 2.1 x 1(T 2 2.1 x 10~ 2 2.0 x 1(T 2 

87 Zr 2.1 x 10~ 2 9.1 x 10- 3 5.1 x 10" 3 3.8 x 10" 3 3.0 x 10" 3 

88 Zr 5.9 x 10" 3 5.5 x 10" 3 5.0 x 10~ 3 5.1 x 10~ 3 5.3 x 10" 3 

89 Nb 1.8 x 10~ 2 2.0 x 10~ 2 1.9 x 10~ 2 1.8 x 10~ 2 1.8 x 10~ 2 

90 Mo 1.2 x 10~ 2 1.1 x 10~ 2 9.9 x 10" 3 9.6 x 10~ 3 9.2 x 10" 3 

91 Nb 8.2 x 10~ 3 9.4 x 10~ 3 1.0 x KT 2 1.0 x 10" 2 9.9 x 10" 3 

92 Mo 7.6 x 10" 3 4.1 x 10 -3 2.3 x 10~ 3 1.8 x 10" 3 1.4 x 10~ 3 

93 Tc 9.5 x 10~ 3 8.8 x 10~ 3 8.1 x 10" 3 7.6 x 10" 3 7.4 x 10~ 3 

94 Tc 2.6 x 10~ 2 4.6 x 10" 2 4.6 x 10~ 2 4.6 x 10~ 2 4.6 x 10~ 2 

95 Ru 3.1 x 10~ 2 2.4 x 10~ 2 1.4 x 10~ 2 9.3 x 10" 3 7.6 x 10" 3 

96 Ru 8.1 x 10" 3 8.7 x 10~ 3 7.4 x 10~ 3 6.6 x 10~ 3 6.2 x 10" 3 

97 Ru 1.0 x 10~ 2 1.2 x 10~ 2 1.0 x 10~ 2 9.5 x 10~ 3 9.1 x 10" 3 

1.2 x 10~ 2 1.4 x 10~ 2 1.3 x KT 2 1.3 x 10~ 2 1.2 x 10~ 2 

1.3 x 10~ 2 2.1 x 10~ 2 2.1 x 10~ 2 2.0 x 10~ 2 1.9 x 10~ 2 

1.6 x 10~ 2 9.2 x 10~ 3 8.7 x 10~ 3 8.0 x 10~ 3 8.2 x 10~ 3 

1.7 xlO- 2 1.8 x 10~ 2 1.8 x 10~ 2 1.8 x 10~ 2 1.7 x 10~ 2 
1.6 x 10~ 2 1.8 x 10~ 2 1.9 x 10~ 2 1.9 x 10~ 2 1.8 



98 Ru 

"Rh 1 Q v m-2 o i v in-2 o 1 v m-2 o n v m-2 1 o v in-2 

ioo pd 

ioi pd 

102p d 
103 Ag 
104p d 

105 Ag 
106 Cd 

107 Cd 



1.6 x 10~ 2 1.8 x 10~ 2 1.9 x 10~ 2 1.9 x 10~ 2 1.8 x 10~ 2 

2.6 x 10~ 2 3.2 x 10~ 2 3.2 x 10~ 2 3.3 x 10~ 2 3.2 x 10~ 2 

5.3 x 10~ 2 7.2 x 10~ 2 7.4 x 10~ 2 7.5 x 10~ 2 7.5 x 10~ 2 

7.5 x 10~ 2 1.3 x 10- 1 1.4 xlO- 1 1.4 x lO" 1 1.4 x 10" 1 

1.2 x 10~ 2 9.8 x 10" 3 5.3 x KT 3 3.2 x 10" 3 2.5 x 10" 3 

9.5 x 10~ 3 6.4 x 10~ 3 3.1 x 10~ 3 1.9 x 10~ 3 1.5 x 10~ 3 
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Fig. 1. — Main nuclear activity at the innermost envelope shell for model 1 (Mns = 1.4 M , 
M acc = 1.75 x KT 9 M yr -1 , Z = 0.02), at the early stages of accretion ( T base = 9.9 x 10 7 
K). Upper panel: mass fractions of the most abundant species (X > 10~ 5 ); Lower panel: 
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Fig. 2.— Same as Fig. HI but for T base = 2.1 x 10 8 K. 
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Fig. 3.— Same as Fig. HI but for T base = 4 x 10 8 K. 
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Fig. 4.— Same as Fig. EEJ but for T base = 7 x 10 8 K. 
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Fig. 5.— Same as Fig. [fl but for T base = 9 x 10 8 K. 
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Fig. 6.— Same as Fig. EEJ but for T base = 10 9 K. 
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Fig. 7. — Time evolution of density (upper left panel), temperature (upper right), pressure 
(lower left), and nuclear energy generation rate (lower right), at the innermost envelope 
shell for model 1 (M NS = 1.4 M , M acc = 1.75 x 1(T 9 M yr" 1 , Z = 0.02), along the first 
bursting episode. The origin of the time coordinate is arbitrarily chosen as 21,150 s, for 
which T base ~ 3 x 10 8 K. 
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Fig. 8. — Same as Fig. but for the overall neutron star luminosity (left panel), and 
envelope size (right panel), as measured from the core-envelope interface. 
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Fig. 9. — Same as Fig. [TJ but for the time when temperature at the envelope base reaches 
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Fig. 10.— Same as Fig. OQ but for T base = 9 x 10 8 K. 
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Fig. 11.— Same as Fig. OQ but for T base = 8 x 10 8 K. 
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Fig. 12.— Same as Fig. OQ but for T base = 4.3 x 10 8 K. 
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Fig. 13.— Same as fig. EQ but for T 6ase = 2 x 10 8 K. 
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Fig. 14. — Same as Fig. (TJ but for the time when temperature at the envelope base achieves 
a minimum value of T m « n = 1.67 x 10 8 K. 
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Fig. 15. — Profiles of density (upper left panel), temperature (upper right), nuclear energy 
generation rate (lower left), and envelope size (measured from the core-envelope interface; 
lower right panel), for model 1 (M NS = 1.4 M , M acc = 1.75 x 1(T 9 M yr -1 , Z = 0.02), 
along the first bursting episode. Labels indicate different moments during the TNR, for 
which T base reaches a value of: (1) 9.9 x 10 7 K, (2) 2.1 x 10 8 K, (3) 4 x 10 8 K, (4) 6 x 10 8 
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Fig. 16. — Left panel: same as Fig. [151 but for pressure. Right panel: Time evolution of the 
total pressure at the outermost envelope shell. 
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Fig. 17. — Left panels: mass fractions of the ten most abundant, stable (or r > 1 hr) 
isotopes for model 1, at the end of the first (upper panel) and fourth burst (lower panel), 
respectively. Right panels: same as left panels, but for overproduction factors relative to 
solar (for / > 10~ 5 ). 
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Fig. 18. — Left panel: Schematic representation of the convective regions in the accreted 
envelope for the four bursts computed in model 1. To simultaneously display the convective 
stages of the TNRs as compared with the overall duration of the bursts, a discontinuous 
time axis is used. Right panel: Temperature profiles, showing the location of the ignition 
point (~ 5.6 m above the core-envelope interface), along the fourth bursting episode, for 
model 1 (M NS = 1.4 M , M acc = 1.75 x 1(T 9 M yr _1 , Z = 0.02). Labels indicate different 
moments during the TNR, for which the temperature at the ignition shell reaches a value of: 
(1) 3.4 x 10 8 K, (2) 4.5 x 10 8 K, (3) 7 x 10 8 K, (4) 9.2 x 10 8 K, (5) 9.7 x 10 8 K, (6) 1.04 x 10 9 
K, (7) 1.06 x 10 9 K, (8) 8.9 x 10 8 K, (9) 7.4 x 10 8 K, and (10) 5.3 x 10 8 K. 
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Fig. 19. — Mean post-burst composition in the envelope for each of the bursting episodes 
computed for model 1 (left) and model 2 (right). 
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Fig. 20. — Light curves corresponding to the second (upper left panel), third (upper right), 
and fourth bursts (lower left), and for the overall computed time (lower right), for model 1. 
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Fig. 21. — Upper panel: main nuclear activity at the ignition shell (~5.6 m above the core- 
envelope interface), when temperature reaches a peak value of T peak = 1.16 x 10 9 K, during 
the 4 th burst computed for model 1. Lower panel: main reaction fluxes (F > 10~ 9 reactions 
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Fig. 22. — Main nuclear activity at the innermost envelope shell for model 3 (Mns = 1.4 
M Q , M acc = 1.75 x 10 -9 M yr _1 , Z = 10 -3 ), at the time when temperature at the envelope 
base reaches T 6ase = 10 9 K. Upper panel: mass fractions of the most abundant species (X 
> 10~ 5 ); Lower panel: main reaction fluxes (F > 10~ n reactions s _1 cm" 3 ) responsible for 
the nuclear activity in the A=65-100 mass region -except for equilibrium (p, 7)-(7, p) pairs. 
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Fig. 23.— Same as Fig. E2J but for T base = 1.3 x 10 9 K. 
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Fig. 24. — Same as Fig. [22j but for the time when temperature at the envelope base reaches 
a peak value of T pea k = 1.4 x 10 9 K. 
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Fig. 25.— Same as Fig. [221 but for T 6ase = 1.3 x 10 9 K. 
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Fig. 26.— Same as Fig. [221 but for T 6ase = 1.2 x 10 9 K. 
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Fig. 27.— Same as Fig. M, but for T base = 10 9 K. 
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Fig. 28.— Same as Fig. [221 but for T base = 7.6 x 10 8 K. 
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Fig. 29. — Same as Fig. [22j but for the time when temperature at the envelope base achieves 
a minimum value of T m ,„ = 2 x 10 8 K. 
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Fig. 30. — Upper left panel: mass fractions of the ten most abundant stable (or r > 1 
hr) isotopes, at the end of the first burst, for model 3. Lower left panel: mean post-burst 
composition in the envelope for each of the bursting episodes computed for model 3. Right 
panels: overproduction factors relative to solar (/ > 10~ 5 ), at the end of the first (upper 
panel) and third bursts (lower panel), for model 3. 
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Fig. 33. — Upper panel: main nuclear activity at the ignition shell when temperature reaches 
a peak value of T peak = 1.23 x 10 9 K, during the 4 th burst computed for model 3. Lower 
panel: same as in the upper panel, but at the end of the burst. 
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episode of model 3. 



